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INTRODUCTION 


The objective of this research participation In the MAGSAT program 
has been to provide experimental and analytical data on magnetic 
mineralogy as an aid to the Interpretation of magnetic anomaly maps. 

The approach has been an Integrated program, ranging from the 
chemistry of materials from 100 or more km depth within the Earth, to an 
examination of the MAGSAT anomaly maps at about 400 km above the Earth's 
surface. Within this framework we have concentrated on providing a 
detailed picture of the pertinent mlneralogical and magnetic relation- 
ships for the region of West Africa, but generalizing the Information 
where justifiable. 

Specifically, our effort has been directed toward 1) examining in 
detail the geochemistry, mineralogy, magnetic properties, and phase 
relations of magnetic oxides and metal alloys in rocks demonstrated to 
have originated in the lower crust or upper mantle, 2) examining in 
general terms the assumption that these rocks, as exposed and collected 
at the surface, portray the nature of their source regions, and 3) exam- 
ining in detail the region*.! geology, tectonics, gravity field and the 
MAGSAT anomaly maps for West Africa. 

The following six-part section outlines the results of the program. 
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FINAL REPORT OF RESEARCH RESULTS 

1 ) OXIOATION PROFILING OF THE CRUST AND MANTLE 
Significance 

Iron In some form Is commonly believed to be dominantly responsible for 
magnetization resident In the outer layers of the Earth, but magnetic minerals 
containing Iron are stable only at specific levels of oxidation so that the 
magnetization of the Earth's lithosphere depends not only on composition, 
temperature and total pressure but also on the partial pressure of oxygen 
fugaclty (fOg). Therefore, an understanding of the oxidation state of the 
lithosphere Is crucial to any mlneraloglcal Interpretation of the origin of 
magnetic- anomalies measured by MAGSAT. 

Resul ts 

A graphical approach to the expected ranges of oxidation as a function 
of temperature and pressure is taken in Fig. 1, where the enstatite + mag- 
nesite « olivine + graphite (EMOG) oxygen buffer is based on Eggler and Baker 
(1982). We have included in our plot the basalt-eclogite transition, the 
shield and oceanic geotherms, and have also calculated and projected the buffer 
curves fayallte ■ magnetite + quartz (FMQ) and wdstlte - magnetite (WM) that 
are more commonly applicable to 1 atm. conditions (i.e. at the surface of the 
earth). Becau.se all of the curves are relatively steep, the relative changes 
in oxidation state with depth do not change dramatically along any one buffer. 
However, the progressive changes that do apply are along the geothermal grad- 
ients; a substantial change in oxidation state toward a more oxidized condi- 
tion of each buffer apparently takes place with increasing depth. 

Fig. 2 illustrates the results of an analytical approach to the question 
of oxidation state at depth, where we have compiled and updated (from Haggerty, 
1976) temperature and jFOg data for basalts; these are shown to have equilibrat- 
ed along or below FMQ. On comparing these data with the oxidation profiling 
shown In Fig. 1 we find that If the depth estimates at which basalts are con- 
sidered to originate (i.e. 20-100 km) are used on Fig. 1, then the levels of 
oxidation that are preserved in basalts extruded at the surface are comparable 
to the oxidation levels calculated. This Is encouraging because it suggests, 
for basalts at the least, that the analytical techniques employed in generat- 
ing the data compiled In Fig. 2 provide reliable probes to the oxidation state 
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Figure 1. Calculated buffers In pressure-temperature coordinates for the 
reactions enstatlte + magnesite * olivine + graphite (CMOS), wiistlte - mag- 
netite (WM) and fayallte * olivine + quart 2 (FMO) illustrating the expected 
ranges of oxidation state for upper mantle conditions. Isobarlc oxygen 
fugacltles for each buffer are plotted within the figure In units of -log 

ffr 








Figure 2a. Sunmary of temperature and fpg data for terrestrial oceanic 
island and ocean floor basalts. Data fall within the hachured envelope; 
reference buffer curves are fayalite ■ magnetite + quartz (FMO) and mag- 
netite - wdstite (MW). 











Figure 2b. Summary of temperature and fO^ data for terrestrial continental 
basalts. Details as In Figure 2a. 
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of the rocks' source regions. 

In further application of the particular technique of evaluating T and 
f0 2 by analyzing the exchange equilibria between Fe and T1 and between Fe 2+ 
and Fe 3+ In co-equIlibrated llmenlte and spinel, we examined llmenltes from 
kimberlites In west Africa. We then compared the results to our compilation 
of T and f0 2 data on other rocks, derived by the same and other methods. Full 
details of the analyses and the compilation are given by Haggerty and Tompkins 
(1983; Appendix 1), and the principle result of this oxidation profiling Is 
shown In Fig. 3. 

The result Is a series of models presenting prevalent oxygen fugacltles 
and rock types for typical cross-sections through the crust and mantle accord- 
ing to tectonic regime. The Earth appears to be stratified with respect to 
redox, such that the surface and near-surface are characterized by magnetite - 
hematite (MH) to FMQ conditions; below this Is a geochemically depleted layer 
within the lithosphere which Is of variable thickness and depth but which Is 
at or below wiistl te-magnetl tr< (VIM), and below this Is a geochemically fertile 
asthenosphere which Is relatively oxidized and characterized by FMQ to WM 
conditions at the shallowest depths, ranging down to wdstlte-iron (WI) and to 
iron at the mantle-core interface. 

Summary 

We have shown that computed and analytical profiling of the oxidation 
state of the crust and mantle are compatible, and that these regions of the 
Earth are stratified in redox potential, the upper crust and the geochemically 
fertile asthenosphere being relatively oxidized and sandwiching a depleted 
and reduced lithosphere. The oxidation state models of Fig. 3 now permit 
meaningful estimates of the potential stability of Iron-bearing magnetic min- 
erals In any selected regime of the lithosphere, and therefore constrain min- 
eralogical models of MAGSAT-measured magnetic anomalies. 

2) SERPENTINITES AS MAGNETIC ANOMALY SOURCES 
Significance 

Whereas it is generally accepted that serpentinltes are possible con- 
tributors to crustal magnetization and that pure magnetite is overwhelmingly 
responsible for crustal magnetization, the magnetic phase equilibria of ser- 
pentinltes and the contribution of metal alloys to serpentinltes' magnetiza- 
tion and the long-range geophysical signature of serpentinltes are poorly 


Figure 3. A diagrammatic cross-section through the crust and upper mantle 
for a variety of tectonic settings, Illustrating the distribution of rock 
types and profiling of oxygen fugaclty. The highest oxidizing conditions are 
represented by the magnetite - hematite (MH) buffer; with Increasing re- 
duction the buffer equilibria are NNO (N1-N10), FMQ (fayallte - magnetite 
- quartz), WM (wdstlte - magnetite) and IW (iron - wflst1te]_. From Haggerty 
and Tompkins (1983). 
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established, to date. 

Results 

We have undertaken. In a companion study, a set of experiments related to 
the stability of Fe-Nl metal alloys as a function of temperature and oxidation 
state. A summary of these data (Fig. 4) Illustrates the regimes of mag- 
netite and metal alloy stability and demonstrates that these phases are po- 
tential sources of magnetism at depth, subject to geothermal gradient. More- 
over, the data are entirely consistent with those Illustrated In Figs. 1-3. 

We have also conducted a very detailed analytical study of serpentlnlzed 
perldotltes from the Josephine Creek area, Oregon, U.S.A. A catalog of electron 
microbeam analyses, rock descriptions, densities, and degrees of serpent- 
inization, analytically establishing that NigFe (awarulte} and CoFe (wairaulte) 
can coexist with magnetic oxides at oxidation levels compatible with ser- 
pentinlzatlon In the upper mantle. Is presented in Progress Report 10/26/1982. 
Microbeam analyses of metal alloys are summarized in Fig. 5., which Illus- 
trates the dominant magnetic alloy compositions and the Curie temperatures 
to be expected. 

The majority of the analysed awarultes from In situ Josephine Creek per- 
ldotltes have Nl/Fe ratios close to 3:1, Indication Curie temperatures of 
580°c + 5°c, but the average Co content is about 1 .25 wt % and the Curie temp- 
erature Increases by about 5°c per wt % Co; additionally, a small but sigi- 
nlf leant percentage of analyses in the range 60-70 wt % N1 Indicate Curie 
temperatures In excess of 590°c even disregarding the affect of Co In awarulte. 
The analyses of walrauites show Co/Fe ratios generally less than the ideal of 
1:1, and mostly Indicate Curie temperatures In the range 940-980°c; Ni con--, 
tents ranging up to 3 wt % may decrease the Curie temperatures by 20-25°e, 

Because remote magnetic measurements alone are equivocal as regards the 
source lithology we attempted to constrain more closely the expected geophy- 
sical signature of serpentinltes by measuring both the magnetic susceptibility 
and density of the Josephine Creek serpentinltes. The result Is summarized in 
Table 1 and Fig. 6. 

Fig. 6 shows the results of our analyses in the context of a graphical 
compilation of data from the literature outlining the relationship between 
magnetic susceptibility and degree of serpentinizatlon for a variety of ultra- 
mafic rocks. From Table 1 and Fig. 6 we draw the following conclusions: 1) 

Magnetic susceptibility and density of Josephine Creek perldotltes Increase 






L 


Figure 4. Log f0 2 vs temperature diagram for the Fe-‘Ji-0 system (from 
MacMahon, 1984). Note especially the T-fp 2 region in which magnetite 
(so 7 id solution) plus Fe-Ni alloy co-exist. 






‘Figure 5a. Frequency- comp os i tion histogram of in situ awaruites from 
Josephine Creek peridotites. The percentages of analyses are relative 
to the number of analyses per hand sample. Composition is presented as 
wt X Ni in Ni-Fe alloy, calculated from the analyses but assuming zero 
per cent impurities. Co is always present, averaging 1.25 wt X and 
ranging up to 2.5 wt X, increasing the effective Curie temperatures. 

S and Cu may also be present but are almost always less than 0.1 and 
1 .0 wt % respectively. 

Figure 5b. Frequency-composition histogram of in situ and nugget 
wairauites from Josephine Creek peridotites. Calculated similarly to 
Figure la. Ni is present in the range 0.6 to 3.2 wt X, slightly de- 
ceasing the Curie temperatures , but other impurities (S and Cu) are 
less than 1 wt % when combined. 
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TABLE 1. JOSEPHINE CREEK PERIDOTITES 
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Hie range of values tor a given sample is due primarily to uncertainty m thi 
fonn of the self-demagnetizing field; samples were measured on a Bison susc- 
eptibility bridge; k c T “4H(k ). 

CgS 

Magnetite is' present in all samples 




Figure 6. Serpentini 2 ation trends in susceptibility-density. Values 
for the Josephine Creek peridotites are in Table 1. Percentages shown 
on the figure are the degrees of serpentinization, with good general 
agreement between the data sets. The theoretical curves calculated by 
Henkel (1976) for the serpentinization of olivines assumed a 100 % 
efficiency in converting Fe released from olivine into magnetite. The 
curve for orthopyroxenes is calculated from molar susceptibility values 
given by Nagata (1961) and handbook values for densities using the 
relation k vo j. C g S * (densi ty) .(molar susceptibili ty)/ (molecular weight). 
The opx and ol curves effecti vely outl ine the allowed region for un- 
serpentinized bimineralic haryhurgi tes . The trends from Henkel (1976) 
and lienert and Wasilewski, 1979) place likely upper and lower limits 
on magnetic susceptibility for a range of serpentini zed ultramafics. 

The curve for Red Mountain peridotite is calculated from mathematical 
and graphical relationships in Saad (1969). Metal alloys begin to 
appear in these Josephine Creek peridotites at serpentinization values 

exceeding 25 %, and are notably present, along with pyrrhotite or 
chromite, in the samples which depart from the general concave-up trend. 



SERPENTINIZATON TRENDS IN DENSITY - SUSCEPTIBILiTY 
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and decrease respectively with the degree of sepentlnlzatlon, conforming to 

general trends culled from the literature. 2) Susceptibilities of the ser- 

-3 .2 

pentlnlzed harzburgltes approach or are In excess of 10 cgs (10 S.I.) at 

about 75 S serpentlnlzatlon and a density of about 2.6 gm/cm^. 3} Josephine 
Creek perldotlte samples which depart from the general trend have significant 
amounts of other magnetic m1nerals« e.g. metal alloy (awarulte and walraulte), 
pyrrhotlte or chromite. 

Stannary 

We have Illustrated the T and fOg stability of co-existing Fe-Ni alloy 
and magnetite* and find the stability field to be consistent with models pre- 
sented In the previous section. We have described the prevalent metal alloy 
compositions present in the Josephine Creek serpentinites, and have noted Ni- 
Fe compositions not conforming to the awarulte stoichiometry. We have defined 
the susceptibility-density field for serpentinlzed ultra-mafic rocks, and 
have shown that the geophysical signature of such rocks is a strong magneti- 
zation associated with a low density. 

3) MAGNETIC SPINELS AND RELATED PHASES FROM MANTLE-DERIVED KIMBERLITES 
Significance 

It Is frequently assumed that the Earth's mantle Is non-magnetlc and 
that lithospheric magnetization resides entirely in the crust* above the 
seismic Moho. However, detailed examinations of mantle-derived materials 
specifically with a view to the presence and origin of opaque phases and 
their magnetic significance are few and far between. 

Results 

A detailed petrographic and electron microbeam study of intrusive kim- 
berlite dike materials from the Koldu complex, Sierra Leone has revealed 
high concentrations of spinels and magnetite. The materials are fine grained 
and In the case of the dikes Intrusion depths are placed at 2-3 km. From 
the silicate phase equilibria and the association of diamond-bearing nodules 
the kimberlites demonstrably originated at depths In excess of 100 km. The 
depths of equilibration of Fe-Ti oxides In these rocks remains to be precisely 
determined but must have been consistent with the MW buffer (see Fig. 5 in 
Appendix 1). The equilibration trends of magnesian ilmenites clearly show 
that ferrimagnetic phases approaching Fe^ In composition result from sub- 
. solidus reactions. The analytical data are deferred to Appendix 2 and are 
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sunmarlzed In Fig. 7, where the significant trend to be noted Is from Cr-rlch 
spinels at the base of the prism toward FejO^ at the apex. 

In the broader context of Iron-related phase equilibria and the Import- 
ance of oxygen fugaclty, we also report the existence of a phase which appears 
to be a new mineral, designated mineral "H" (Tompkins, 1983). Compositions 
are detailed In Appendix 2 and sunmarlzed In Fig. 8. The mineral is con- 
sidered to be analogous to a Mg-A1-T1-Cr-w<lst1te-type mineral. Oxygen fu- 
gacltles for formation of this mineral are estimated to be low, on the order 
of 10* 11 to Itf 12 atms at H00°c, but the magnetic significance of this phase 
remains to be established. X-ray data are required In order to establish 
the crystal structure and the magnetic ordering of the phase, but small grain 
sizes (» 0.05 mm) have so far prevented extraction of discrete grains. 

A further indication of the possible magnetic significance of mineral 
"H" Is drawn from paleomagnetlc experiments on cores from Dike Zone B at 
Koldu (Tompkins 1983). Iron-bearing phases observed In these specimens in- 
clude llmenites, spinels, pyrrhotlte, mineral "H" and Iron metal. On a 
Zljderveld plot (Fig. 9a) there are two components of magnetization, over- 
lapping In their coercive force specta. The early removal of a soft com- 
ponent is attributed to Iron metal and mineral "H", the latter expected to 
have a coercive force Intermediate between that of Iron metal and magnetite. 
The decay of magnetic Intensity (Fig. 9c) of core 403 is also consistent 
with the early removal of a low coerclvlty component. 

Suitmary 

We have shown 1) that klmberlltlc llmenites are capable of breaking down 
to form ferrlmagnetic spinels close to magnetite In composition, and that this 
occurs at depths between *100 km and the Earth's surface, and 2) that de- 
tailed petrographic, electron microbeam and paleomagnetlc studies implicate 
highly reduced Iron minerals as significant in kimberlite magnetism. 

4) MAGNETITE AND MAGNETIC SUSCEPTIBILITY RELATIONSHIPS 
Significance 

In order to either Interpret or model magnetic anomalies on a mineral - 
oglcal basis and as due to induced magnetization it is necessary to quantify 
the relationship between mangetic mineral content and magnetic susceptibility. 
Furthermore, if metal alloys are to be invoked as partially responsible for 
magnetization,- it Is vital to first understand the role to be played by mag- 


Figure 7. Compositional spinel trends for Koldu, summarized from Appendix 2, 
plotted on the oxidized spinel prism and compared with spinel trends from 
southern Africa and Canadian shield. Mt * magnetite; Mag.fer. * nagneslo- 
ferrlte; Her ■ hercynlte; Sp ■ spinel; Chr « chromite; Plc.chr. » plcro- 
chromlte. Klmberlltlc trend 1 Is from Mitchell and Clarke (1967), and trends 
2:and'3 are •from Haggerty (1976);' Note especially that the trend for Koldu 
spinels Is from the base of the prism toward compositions close to Fe 3 0 4 at 
the apex. 










Figure 8. Mineral "H" Is plotted as a function of atomic proportions in 
terms of T1, and R^ + . The field of wlistite (Fe-|_ x 0) is outlined In the 
hatched area, and the spinel Is drawn connecting endmember ulvdspinel ("Usp") 
to magnetite (Mt) or chromite (Chr). Stoichiometry of mineral "H" Is re- 
presented as MgOy. See Appendix 2 for analyses of mineral "H". 



Figure 9. Alternating field demagnetization data to 50 MT for Dike Zone B 
(samples 401, 402, 403). (A) Zljderveld plot, with the area around the 

origin enlarged In A. Open circles ■ vertical component of magnetization 
(E-W versus U-D) ; solid circles ■ horizontal component of magnetization 
(E-W versus N-S); W *west; E ■ east; N * north; S * south; U « up; D * 
down. (B) demagnetization decay curve. (C) demagnetization decay curve 
normalized to NRM intensity. 
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netlte. In addition, as noted above. Iron oxides and alloys are subject to 
redox conditions. This section evaluates the relationship between the oxida- 
tion state and amount of magnetite present In a rock and the resultant mag- 
netic susceptibility. 

Results 

Fig. 10 presents the results of a literature search and compilation of 
relationships between volume 5 magnetite versus room temperature susceptibil- 
ity. Although there Is a large degree of scatter within each data set and the 
relation Is plotted In log-log space, the good overall agreement between diff- 
erent sets of natural samples and between the empirical and theoretical re- 
lations provides a guide for magnetic anomaly Interpretation based oh In- 
duced magnetization. On this basis, Induced magnetizations In the range 1-10 
A/m are generated in the Earth's field by 1-10 volume per cent magnetite, 
ranges normally attained only by basaltic rocks and iron-formations. 

However, this apparently useful and consistent relationship Is mislead- 
ing In two respects, as illustrated In F1g.il which shows In more detail the 
data compiled by Dunlop (1974), simplified in Fig.io as a straight line. 

Fig. li demonstrates 1} that the room temperature susceptibility of a given 
rock type may vary over a wide range, and 2) that for basalts, and by im- 
plication other rocks, the susceptibility depends not only on the amount of 
opaque oxide minerals generally identified as magnetite, but also depends 
critically upon the oxidation state of these minerals. This Is dearly ill- 
ustrated In the top right of the figure, where, for a given amount of "mag- 
netite", the susceptibility varies according to oxidation state, such that the 
more highly oxidized samples have lower magnetic susceptibilities. 

Summary 

We have confirmed that there exists a first order correspondence between 
magnetite content and magnetic susceptibility and that this may be loosely 

correlated with lithology. Induced magnetizations of 1-10 A.nf 1 require 
1-10 vol % Fe 3 0 4 at room temperature. Most significantly, in view of pre- 
vious discussion as to the redox state of the lithosphere, we have pointed 
out that oxidation state is a critical parameter affecting susceptibility, 
by virtue of the redox response of Fe-Ti oxides. 


. — *■*» • * *w«i 


Figure 10. Relationships of apparent magnetic susceptibility with 
volume " magnetite. Empirical relationships are for rocks from the 
Adirondacks (Balsley and Buddington, 1958), Minnesota (Mocney and 
Bleifuss, 1953) and diverse locations (Dunlop, 1974), along with the 
theoretical relationship of Stacey and Bannerjee (1974), a two-fold 
thermally enhanced susceptibility (Dunlop, 1974), and typical average 
susceptibility ranges (Dobrin, 1976; Nagata, 1961). Named rocks are 
from Minnesota (Mooney and Bleifuss, 1953). Mantle xenolith specific 
susceptibi 1 ties (Wasilewski et al . , 1979), recalculated to volume 
susceptibilities on the basis of average densities from Daly et al. 
(1966) and Robie et al , (1966), are shown with horizontal bars. The 
susceptibility ranges for induced magnetizations of 1 A/m and 10 A/m 
over the range of intensity of the Earth's magnetic field are also 
illustrated. Note that the left and right ordinate scales are equi- 
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VOLUME % MAGNETITE 


Figure 11. Apparent volume magnetic susceptibility versus volume % 
magnetite, plotted from a compilation of data by Ounlop (1974). 

Roman numerals attached to the basalt data points are the oxidation 
indices of Wilson and Watkins (1967), wherein Class I is homogenous 
titanomagnetite, and progressive oxidation is toward Class V pseudo- 
morphing of titanomagnetite by pseudobrookite and hematite or by 
spinel and rutile and hematite. 
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5) MAGSAT VECTOR ANOMALIES RELATED TO THE WEST AFRICAN AND AMAZONIAN CRATUNS 
Significance 

Most Interpretations of the MAGSAT anomaly maps to date have employed 
the total field anomalies. Much of our analytical data on magnetic mineralogy 
and oxidation profiling stems from rocks collected within west Africa. Ue 
have examined the MAGSAT vector anomaly maps for the regions of west Africa 
and northern South Amer1ca t alone and In relation to the regional geology 
and geophysics, drawing upon our field experience and our mlneraloglcal data. 
Results 

The principal result is our manuscript "MAGSAT vector anomalies related 
to the west African and Amazonian cratons", reproduced as Appendix 3. 

Summary 

Following Is the abstact form our manuscript: 

ABSTRACT 

MAGSAT (MGST 4/81) anomalies are similar over the West 
African and Amazonian cratons., Differences between total field 
(aB) and N^S horizontal vector (aX) anomalies over both cratons 
are largely those that would be expected for anomalies due to 
Induced magnetization. At the present low magnetic inclinations, 
the negative AX anomalies which are spatially associated with 
the cratonic shields Imply normally-directed positive magneti- 
zation-contrast sources. Archean nucliee of the shields appear 
to be more strongly magnetized than Proterozoic segments, but 
both stand out from adjacent basins and ford belts. Under spe- 
cified conditions, the uncertainty inherent in the poorly de- 
fined MAGSAT anomaly profiles taken along opposite siedes of the 
continental rejoin In a pre-drift reconstruction. Synthesis model 
anomalies generated form lithological and gravity data roproduce 

reasonably well the cratonic negative MAGSAT anomalies. The 
dominant contribution to MAGSAT aX anomalies over the West 
African and Amazonian cratons results from magnetization con- 
trasts within the crust. These are modulated by the effects of 
structural contrasts, intrusions and metamorphism within the 
lithosphere, and are localized by the preferred distribution, 
of iron deposits into the more ancient crustal blocks. 
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6) NATIVE IRON AND THE OXIDATION STATE OF THE CONTINENTAL LOWER CRUST 
Significance 

Although we have developed models for the gross oxidation profile of the 
upper mantle* the oxidation state of the continental lower crust has not been 
firmly established, tb-date: This knowledge has remained a vital missing 
link, given the dependence of magnetic mineralogy upon oxidation state. An 
integral part of our research related to MAGSAT has been an effort to Iden- 
tify materials of unquestioned lower crustal origin which could shed light 
on this problem. This section describes the results of this search. 

Results 

The principal result Is our recognition of native Iron In lower crustal 
granulites recovered from a Kimberlite pipe In west Africa; the phase 
equilibria and magnetic mineralogy strongly suggest that the lower crust Is 
In a highly reduced state and that Iron metal may contribute to MAGSAT anom- 
alies. Full details of this research are In our manuscript "Native iron In 
the continental lower crust : petrological and geophysical Implications", 
submitted to Science (August, 1984), see Appendix 4. The abstract follows : 

Abstract. Lower crustal granulite xanoliths recovered from a kimberlite pipe 
in west Africa contain native iron (Fe°) as a decomposition product of both 
garnet and ilmenite. Magnetic measurements show that less than 0.1 volume 
per cent of iron metal is present. Geothermometry and oxygen geobarcmetry 
indicates that the oxide and metal phases equilibrated between iron-wustite 
and magnetite-wustite buffers which may represent the oxidation state of the 
continental lower crust, and the depleted lithospheric upper mantle. Ferro- 
magnetic Fe° could be stable to a depth of «95 km, and should, be considered 

in the interpretation of long wavelength static magnetic anomalies. 

Summary 

We have shown that native Iron Is present In some lower crustal gran- 
ulites, and that the continental lower crust Is probably at low oxidation 
levels compatible with metal and metal alloy stability fields. The impli- 
cation for MAGSAT is that the interpretation of anomalies needs not, in fact 
should not, be constrained as due to magnetite. The magnetic outer shell of 
the Earth may be as thick as »100 km in places and Iron metal, and perhaps 
other metal alloys, should be considered in any interpretation of MAGSAT 
magnetic anomalies. 
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SUMMARY 

We have provided extensive data on magnetic minerals and their host 
rocks, particularly for serpentlnlzed ultramaflc rocks, lower crustal meta- 
morphlc rocks and klmberlltlc suites. We have established oxidation pro* 
files for the crust and mantle, and we have sunmarlzed the results of many 
hundreds of electron microbeam analyses describing the chemistry and phase 
relations of Fe-TI oxides In mantle rocks, and of oxides and alloys In a 
collection of serpentlnlzed perldotltes and a collection of lower crustal 
granulltes. 

Overall, we conclude from our analytical and experimental data that 
magnetization need not be confined to the seismic crust and that magnetic 
oxide and alloy compositions are geochemical 1y stable under lower crustal 
and mantle conditions; although It Is not necessary for modelling purposes 
alone to Invoke a mantle source for any particular MAGSAT anomaly, our 
data demonstrate what the hypothesis that only the crust Is magnetized and 
that only magnetite Is responsible for magnetization Is simply an assumption 
of convenience, and one which Is misleading. 

We have analyzed qualitatively the MAGSAT maps over shield regions In 
west Africa and northern South America, and have demonstrated very similar 
relationships for both regions. Obvious among the results from this work 
Is the observation that large deposits of iron ore are visible In the MAGSAT 
maps. 

In conclusion, It is re-emphasized that magnetic anomaly' Interpretations 
is dependant upon knowledge of the mineralogy, and that the MAGSAT anomaly 
maps demonstrably represent the Integrated effects of magnetic minerals with- 
in the mantle, within the crust, and exposed at the very surface of the 
Earth. The two overriding constraints In our analysis are that temperature 
and oxygen fugaclty ( and hence composition) remain the most significant model - 
limiting parameters In the Interpretation of MAGSAT anomalies. 
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diamonds'*. and the sources of deep-seated magnetic 
anomalies" ". Although the study of mantle nodules to deter- 
mine the redos state of the upper mantle is in its infancy there 
has been considerable debate on whether it is reduced 11,1 * and 
is equivalent to equilibrium conditions defined by iron-wuitite 
(IW), or on whether it is more oxidised 1 *' 1 * and, hence, is more 
appropriately defined by equilibria consistent with magnctitc- 
wQstite (MW) or feyeiite-magne tile-quartz (FMQ). In the 
extreme case at magmatic temperatures (1,200*0, these esti- 
mates differ by approximately five orders of magnitude in 
oxygen fugadty (/<*)• Translated, this difference is similar to 
the differences between the oxidation states of lunar and terres- 
trial basalts, a disparity that is unreasonably large. In this study 
we conclude that the source region is compatible with WM- 
FMQ and poasibie explanations for the diversity of redox states 
reported indude the ages of the mantle being sampled and 
whether that portion of the mantle is mature and depleted, or 
juvenile and enriched in volatilea. 

Three methods have been used to estimate the redox condi- 
tion of the upper mantle, and all rely on the variations of iron: 
oxygen ratios in minerals and on the magnitude of Fc’VFc** 
concentrations. Thermodynamic techniques rely on theoretical 
solution models for coexisting minerals ‘‘■* 7 . In the intrinsic fa, 
method, furnace gas mixtures (CO or COj+Hj) are adjusted 
as a function of temperature (T) to match the redox state of 
the sample'^ '*. In the third technique, exchange equilibria 
between Fe and Ti, and between Fe 1 * and Fe J *. define unique 
solutions for T and /©, in co-cquilibratcd members of the 
ilmenite and spinel mineral series’*. The present study used 
the last of these methods and it has also been used to 
estimate the redox states of planetary basalta’*. A critical factor 
in this application it that the mineral pair ia demonstrated to 
have co-equiiibrated, a situation that is difficult to prove because 
ilmenite end spinel members are typically discrete constituents 
in igneous and metamorphic rocks. In the assemblages 
described here, however, the two components are in the grain 
boundary contact, and ionic exchange between ilmenite and 
spinel, in response to external environmental conditions, is 
reassured. 


Samples and analytical techniques 

One aspect of an extensive geochemical programme on the 
characterization of the West African kimberlites and their 
associated upper-mantle derived nodule suites’* 40 , has been a 
detailed investigation of the oxide minerals for exploration 
purposes, magnetic properties and because of their sensitivity 
to modification in oxidizing or reducing conditions. Ilmenite 
occurs in great abundance in the Cretaceous kimberlite fields 
of Liberia, Sierra Leone and Guinea, and ilmenite is part of 
the discrete nodule suite typical of other kimberlite provinces, 
but debated in terms of their precise depths of origin 2 ’. That 
ilmenite nodules in kimberlites originate in the upper mantle 
is indisputable. This is borne out by diamond inclusion 
studies 11 - 33 , by garnet and pyroxene intergrowths 214 *, and 
by high-pressure, high-temperature phase-equilibria 
experiments 2 * 4 *, that all include ilmenite as a component. 

The ilmenites we have examined are typically 1-2.5 cm in 
diameter with some as large as 15 cm. Spherical or oblate 
spheroidal in shape, the ilmenites are free of surface weathering 
but occasionally show some reaction with the enclosing kimber- 
lite, manifested by very thin veneers of late stage perovikitc 
and associated calrite. Examined in reflected light and under 
oil-immersion microscopy, a large number of ilmenites are 
observed to contain parallel, oriented lamellae (20-100 in 
width) of spinet (titanomagnetite) that are crystallographically 
controlled along (0001} rhombohedral ilmenite planes. Among 
these samples are two contrasting types of assemblages. In some 
ilmenites from the Koidu Complex, Sierra Leone, exsolution 
of ilmenite lamellae within the host ilmenite has taken place, 
and this process is followed, or accompanied by subsolidus 
equilibration that has also formed spinel (Fig. la). Rounded 
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Fig. 1 Reflected light photomicrographs taken under oil-immer- 
sion objectives. The width of each plate ia 0.475 mm. The field 
of view shown in a is typical of assemblages observed in ilmenita 
mtgacrysts from Koidu. Sierra Leone. The light grey, host ilmenite 
is the Origins! (W ilmenite that has undergone partial re-equilibra- 
tion to ilmenite (1)+ spinel The bright end highly reflective areas 
are carbonate i ndusions or minor sulphides. In b the «* wk oriented 
lamellae arc spinel in an ilmenite host typical of mcgacrysti from 
Liberie, and the Antoschka kimberlite pipg in Guinea. 

or curvilinear liquid immiscible sulphide (pyrite, pyrrhotite, 
pentlandite, chalcopyrife) and carbonate inclusions are present, 
and the spinels are interpreted to derive from subsotidus reduc- 
tion. In other ilmenites from Koidu, and in all of the ilmenites 
from Liberia and Guinea (Antoschka pipe), the textural rela- 
tions consist simply of parallel sets of spinel lenses elongated 
along {0001} ilmenite planes (Fig. 15). These assemblages also 
result from reequilibration, and it is assumed, based on experi- 
mental studies’* 4 *, that solid state reconstitution his taken 
place in reducing conditions, transforming homogeneous 
ilmenite into a compositionally modified ilmenite + spinel. This 
textural and mineral relation it described as reduction exsol- 
ution. -Although we cannot demonstrate that reduction took 
place in the source region, this appears most probable, and the 
process was deuteric modification with prolonged cooling. We 
can confidently eliminate the imiptive kimberlite event as the 
reduction agent, because of exsolution iuiuut:ricto)ol ilmenite 
and the rarity of external reactions. 

Compositions of coexisting spinel and ilmenites were deter- 
mined by electron microbeam techniques using an ETEC 
Autoprobe, natural and synthetic mineral standards, and con- 
ventional correction procedures 274 *, with standard operating 
conditions (15 kV and 0.03 nA specimen current). Iron(u) con- 
centrations were calculated assuming ideal analyses and mineral 
stoichiometry. Reliable analyses were only possible on the 
coarsest intergrowths, and only those data for which repeated 
and reproducible compositions were obtained are considered. 
Averages presented are a minimum of 10 analyses. 

Uppcr-mantie derived oxides are commonly enriched in 
MgO, MnO, CrjOj, and AljOj, and in the classical application 
of Buddington and Lindsley's 11 Fe-Tt oxide geothermometer 
and oxygen geobarometer, T and fa, estimates were only poss- 
ible in assemblages having negligible contents of elements other 
than Fe-Ti-O. A revised thermodynamic model 3 *, and a 
method of treating these other components has been pro- 
posed 20 , and the 7* and fa, of coexisting oxides in kimberlites 
can now be obtained using the equations: X ,< m » 


After integrating equation < 14) we define 

^ * -0 ~ (15) 

The 8 and 4 component! of the transverse poioidai acceleration 
are respectively' 1 3 V /34 and (iV/A4\fi\nO. The second term 
in the bracket is dominant. The spherical harmonic angular 
dependence of 4’ is also present in f, and in V. Thus the shape 
of the surface distortion ({.) determines the transverse fluid 
velocity from equation (IS). 

Even if there were a non -zero toroidal part of the circulation 
current the 'whole -disk average' of the Doppler velocity would 
praduJc observing it. The gradient and poioidai parts of the 
velocity can be observed, but only the poioidai transverse part 
of a 13-day period can be significant. Choosing the line of tight 
as tltc axis of the spherical coordinate system, only the m -0 
harmonics (Legendre polynomials) can be observed. 

Neglecting limb darkening effects, the symmetry of the whole 
disk avenge imposes other restrictions on the results. Except 
for / • 2, spherical harmonics of even degree are unobservable 
in the whole disk average and /«• 2 yields only the second 
harmonic of the cote rotation frequency. The / ■ 1 degree 
cannot occur in 4'- The l • 3 spherical harmonics yield the 
fundamental and the third harmonic of the core rotation 
frequency. The /«3, n “3 mode (referred to the solar axis) 
is the strongest source of the fundamental period. 

If the above interpretation of the data is correct, the / — 3 
harmonic is probably the source of the circulation current. A 
detailed comparison with the 1966 observations is not useful, 
for the telescope was incapable of detecting odd degree 
spherical harmonics. 

Calculating the 8 component of the velocity for / » 3 from 
equation (15) and integrating over the solar disk gives 


Tj»p(r,) 

2wr*p 


(16) 


for the velocity amplitude where f. • <rf>j(cos 8) and T is the 
core rotational period. Asauming (■3,A/*r-r l «5x JO* cm. 
p(r ( ) ■ 3 u 10\ p • 4 x 10'* and e *■ 10* cm gives 440 cm s*' for 
A. The value of the coefficient a was chosen arbitrarily but it 
is reasonable. 

If the above interpretation it correct there is a surprising 
conclusion. The above value (a) for the odd spherical harmonic 
distortion (f ■ 3) is fairly large, of the same order of magnitude 
•s / « 2 and f *4 distortions 4 . This suggests that if there is a 
strong magnetic Held buried in the solar core, it may be quite 
complex in f orm. 

The centrifugal distortion due to the rotation of the surface 
should also drive an Eddington-Sweet current In this case the 
flow it from the Equator to the pole and the velocity varies as 
tin 2#. The velocity Is limited by turbulent viscosity which is 
not easily estimated. Also this flow may be dominated by other 
effects. 

Durrani and Schrdtcr 11 have suggested that the CUveria tt 
al 13-day period may be a spectroscopic artefact u p da t ed 
with the rotating magnetically active tones on the Sun's surface. 
From the plage maps for the same observational period used 
by Ctaverie et al., they compute a function (here 0(0) which 
is similar in appearance to the ClaVerie function. Their effect 
should exiit, but the magnitude is uncertain. I find that a simple 
sinusoid (Fig. 1) provides a substantially better (robust) fit to 
the velocity data than does the D function. Also, adding the 
D(t) term to the sinusoidal fit does not significantly improve 
the fit Also, the D(0 fits better if it is phase advanced by 2 
days (55*), D -»D(r+ 2). I conclude that the Durrant-Schrdter 
effect should exist but that it is probably weak, and perhaps 
negligible in the Oaverie et al. data, Edmunds and Gough have 
independently considered the same artefact. From a study of 
line formation either side of the potassium line they conclude 
that the magnitude of the effect is approximately right to 
account for the observations 14 . 

This research was supported in part by the NSF. 
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Redox state of Earth’s upper mantle 
from klmberlitic ilmenites 

Stephen E. Haggerty & Linda A. Tompkins 

Department of Geology, University of Massachusetts, Amherst, Massachusetts 01003. USA 


Temperatures and oxygen fugacides are reported on discrete itmenite nodules in kimberlites from West Africa which 
demonstrate that the source region in the upper mantle is moderately oxidised, consistent with other nodule suites in 
kimberlites from southern Africa and the United States. A model is presented for a variety of tectonic settings, proposing 
that the upper mantle is profited in redox potential, oxidized in the fertile asthenosphere but reduced in the depleted 
lithosphere . 


volatiles and the oxidation-reduction environment exerted 
by gas species in the system C-O-H-S are considered to have 
had profound effects on planetary and atmosphere formation' 
and are critical to pctrogenetic models for the generation of 
magmas within the Earth's upper mantle 1 . The compositions 
of volatiles in this source region influences degrees of partial 
melting, the composition of the melt, mass fractionation, liquid 
lines of descent, valency states of elements, and elemental 
partitioning coefficients between liquids tnd solids. Gts samp- 


ling of active subaerial volcanoes 1 * 3 and analyses of trapped 
volatiles in quenched ocean floor glasses have provided a wealth 
of new data on volatiles being generated at present 1 ' 4 * 7 . Samp- 
ling of the upper mantle through nodule suites in alkali basalts 
and kimberlites have received less attention, but the realization 
that these materials provide direct information on the volatile 
budget of the mantle and on the redox state of subcrustal regions 
is rapidly being recognized. These data have far-reaching impli- 
cations ranging from mantle metallogcnesis 4 *. to the origin of 
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(Mf4■0JAI+Fe , *)/(0JFe , *+Mf+04AJ+F« , *). and 
(Ti - 0.5Mn)/(Ti - 0. JMn + 0, JFe r *), where X llm and Afu, arc 
the mole fractions of ilmenite (in ilmcnite solid solution) and 
ulvdspincl (Fe,TiO.) in spinel solid solution, respectively. 
Uncertainties in T and fo, (excluding Cr,Oj) by this method 
are approximately 40-80 *C and 0. 3-1.0 log unit f a „ assuming 
±1% uncertainties in spinel and ilmenite compositions**. We 
have modified the above equations, treating Cr t Oj in a manner 
similar to AltO> and although the precise effect of these oxides 
on T and fo, remains unknown, we estimate that the values 
arc still within the range »100*C and -1 log unit fo,- These 
estimates are corroborated by low CriOj+AJiO, mineral 
pain and by independent experimental data 11 in which we 
consider ilmenite alone in the system FeTiOj(ilmenite)- 
MgTiOitgcikielitei-FejO^haematiU) as a function of fo, at 
constant T (1,300*0. The dau are also consistent with the 
estimates obtained on kimbcriitic Umenitcs from southern 
Africa and the United States, by thermodynamic treatment 17 , 
and using the intrinsic fo, method**. 

Results 

The compositions of 17 iimenite-spinel pairs (live from Koidu 
and six each from Liberia and Antoschka) are expressed as 
major element, combined divalent (FeO+MgO+MnO), 
trivalenf (Fe ,0, + A1,0 j + Cr jO j) , and tetravale nt (TiO j) oxid e$ 
in Fig. 2. Spinels from Koidu lie on or very close to the join 
FejTi0«-R)0« and all remaining spinels, with one exception, 
lie along joins slightly enriched in TiOj, a reflection of their 
MgO concentrations. Partitioning of MgO is preferential into 
ilmenite; hence, all of the ilmenite dau fall between FeTiOj 
and MgTiO). An approximate evaluation of these data, in terms 
of fov (»>y be obtained from the inset in Fig. 2, based on the 
relative positions of coexisting compositions, and the slopes of 
tic-lines between ilmenite and spinel solid solution joins 11 . The 
more enriched Fe** samples (Koidu) are more oxidized (NNO), 
and those with higher Fe** contents (Liberia) are relatively 
more reduced (JVM). 

Absolute T and Am obuined from these dau are illustrated 
in Fig. 3, and four of the standard T-/o> buffers are included 
for reference. These are IW, WM, FMQ, and Mn0-Mn»0 4 . 
Nine of the 17 oxide pairs lie, within experimental limits, along 
the FMQ buffer; live pairs are more oxidized than FMQ; and 
three are more reduced than FMQ. Because the assemblages 
have undergone reduction-equilibration, original single phase 
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Fig. 2 Electron microbeam analyses of coexisting spinel and 
ilmenite compositions as a function of RO (FeO+MgO+MnO), 
RjOj (FejOj + CriOj + AljOj), TiOj. 


ilmcnites must have been stable at Am above those shown for 
each of the coexisting oxide pairs. Also included in Fig. 3 are 
data from the Precambrian Premier kimberlite 11 , ihc 
Cretaceous Oka carbonalite* 4 , the Bcnfontein kimberlite- 
carbonatite sills’* (also Cretaceous in age), experimental data 
on ilmenite from the Monastery (Cretaceous) kimberlite**, 
along with dau from peridoiites enclosed in Tertiary alkali 
basalts’ 1 . The most striking feature of the diagram ii that ail 
data, with the exception of the peridoiites, are more oxidized 
than MW. the peridoiites bracket IW, and these dau, in com- 
mon with other equally reduced values on olivine in basalts”, 
were made using the intrinsic for. Zr-electrolytic cell tech- 
nique**' 17 . Kimbcriitic Umenitcs, however, yield values coin- 
cident with FMQ**. Oxygen fugadty increases with T but also 
increases with pressure, and In Fig. 4 the Ah buffers are correc- 
ted to P <■ 30 kbsr to illustrate this effect The extent to which 
the iimenite-spinel dau are modified at high P it unknown, 
but the sente of the correction it towards higher values of fo,. 
Both the oxide dau and the buffer curves will be synchronously 
dis pl ac ed with increasing P. However, maintaining the oxide 
dau as uncorrected values, because the correction is considered 
to be small 1 *, and assuming 30kbar, the Koidu data approach 
FMQ, the Oka and three of the Liberian compositions are 
equivalent to WM. one it more oxidized and the remainder 
more reduced. The Antoschka data lie -1 log unit below WM. 
Because the precise pressures are not known the relative Am 
are clearly on!.- estimates. 

Kimbcriitic ilr,«enites are dominated by solid solutions among 
three endmember components (FeTiOj-MgTiOj-Fe.Oj) and 
because the oxidation sute of iron in two of these endmembers 
is either exclusively Fe** or Fe’* (that is ilmenite and haematite, 
respectively), the stability of any composition within the 
ilmenite ternary must be fo, dependent. This fo, dependency 
is illustrated in Fig. 5 based on experimental data” at 1,300 *C. 
The West Africa dau are shown together with published dau 
on the compositions of ilmenites from kimberlites, 
worldwide** - * -. Most compositions lie within the range of fo, 
between 10”* and 10"* atm. The WM buffer at 1,300 *C and 
30 kbar, has an Ah * 10~* atm. A more appropriate buffer for 
upper mantle conditions has been proposed by Hggler and 
Baker 1 *, which is defined as EMOG, and it the intersection of 

Temperature < C) 



Fig. 3 Temperature-oxygen fugadty relations for coexisting 
iimenite-spinel pairs from West Africa (Liberia: Antoschka- 
Guinea; and Koidu-Sierra Leone) determined according to the 
methods described in refs 18, 29, 30. Other data shown are for 
the Premier (Prem.) kimberlite”, the Ok* carbonatite 14 , the Ben- 
fontein (Ben) Icimberiite-carbonatite ills 1 ’, the Monastery kirn* 
berlite 1 *, and peridoiites (Per. A and Per. B) in alkali basalts '*■**. 
The solid tines are experimental buffers Mn0-M,0, (ref. 32); 
FMQ (fayalile-tnxgneiiie-quarts)* 1 ; MW (nugnetite-wusiite) and 
IW (iron-wustite) arc from ref. S4. The salient feature to be noted 
it that the peridotile dau are considerably more reduced than 
those from kimberlites and carbon atites. 
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two univirient ructions: (1) enstati te + magnesite * carbon » 
forstcrite+ vapour; with (2) magnetite + carbon **wSitfte + 
vapour. This intartcction is at 29.6 kbar, 10"* atm and 1,230 *C. 
At lower T, EMOG lies between WM and FMQ. but above 
1,230 *C at high pressures, graphite is replaced by diamond 
and EMOD is —1 log unit more reduced than WM. If the 
equilibrated data in Figs 3 and 4 are eslrapolated to 1,300 *C. 
parallel to the FMQ buffer, the range in resultant fo, is 10'** 
to !0~* atm. This corresponds almost exactly to EMOO, WM, 
and the range of compositions in Fig. 3, noting that it is only 
the ilmenite component of the oxide pair that is considered in 
this treatment This estimate in fo% i* consistent with the ther- 
modynamic approach for coequilibrated oxides and silicates 
from a variety of upper-mantle derived nodules 17 -**, and with 
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Fig. 4 Pressure corrected fo, buffers to JO kbar, illustrating the 
effect of P on the relative fo, for the 1 bar data shown in Fig. 3. 
No correction has been applied to the mineral data. 
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Fig. 5 Ilmenite compositions from West Africa (sec legend to 
Fig. 2) plotted as a function of Fe 2 Oj (haematite). FeTiOj 
(ilmenite) and MgTiOj (geikielite). M and a refer to an ilmenite 
Karting composition, and to an experimentally equilibrated com- 
position. respectively from the Monastery kimberlite"*. The gen- 
eralized field for kimbcrli‘ : c ilmenite! worldwide are from a compi- 
lation of dal a cited in refs 23. 38, 39. Oxygen fugacity isobars are 
from ref. 31; WM. wujlite-magnetite 34 and EMOG (see text for 
explanation) it from ref. 16. There is a large sslvus in the ternary 
extending from the FejO.i-MgTiOi sideline towards FeTiOj 
which has been omitted for clarity. It is considered (Hat the 
distribution of kimberlitic ilmenite compositions is influenced hy 
this ternary tolvus and hy the fo, environment of formation 1 ". 


intrinsic /o, data on kimbcriiiic ilmenite". Ths estimate is also 
consistent with the oxygen fugacitics of continental basalts in 
which, *«96% of the data lie between FMQ and WM UI . 

Discussion 

It has been argued* 1 that tha presence of graphite or carbon 
polymorphs, or the occurrences of CH< and CO in ocean floor 
and upper-mantle derived nodules mandates that tha mantle 
be reduced. We consider that these assumptions are erroneous: 

(1) because the abundances of reduced carbon species are 
negligible, amounting in most cases to concentrations of 
■0.5 wt% in media that are dominantly COa+H a O (ref. 2); 

(2) graphite is stable over an enormous range of fo, and 
temperature conditions**; (3) graphite or diamond formation 
is possible by cither oxidation (CHai 5 C+ 2 Hi) or reduction 
(COjSsC+Oi), and diamond is in equilibrium with COj + 
CH«±»2C+2HjO (ref. 10), so that extremely reducing condi- 
tions are not necessarily a prerequisite to carbon precipitation; 

(4) whether carbon forms from a single (such as CO±?C+iOz) 
or a multiple series of reactions (such as COt^CO-f C+iOj), 
in which CO is an intermediate product, typically cannot be 
determined and hence the conditions of /©, are indeterminant; 

(5) high-pressure equilibria in the system C-H-O show that 
CO is never a major C-bearing species in the vapour 10 ; (6) 
methane, however, may be present in mantle conditions 10 
although evidence to the contrary suggests that the equilibria 
are more appropriately defined by CO-COj (refs 16, 17). In a 
broader context wc note that carbon + magnetite are typical of 
chondritlc meteorites and some classes of metamorphic rocks, 
and that these are relatively oxidized. 

Degassing of the upper mantle has taken place with time and 
an inevitable consequence of degassing is that the redox state 
of the mantle must also have changed with the continuous Ibss 
of volatiles. The question is whether the mantle has become 
more reduced or more oxidized, or whether the condition is 
steady state. The latter assumption is that volatile loss through 
volcanism is offset by ocean floor subduction of hydrous 
materials. Present-day basaltic volcanism and the Cretaceous 
kimberlite and carbonatitic suites collectively imply mantle 
redox conditions between FMQ and WM. However, we offer 
two tantalizing and perhaps crucial pieces of data on the redox 
state of the early mantle. These are that Ni-poor metallic Fe 
has been identified as an inclusion in diamond*’ and that the 
ages of diamonds, regardless of the time of kimberlite intrusion 
are Archaean**. Considerably more data are required but we 
propose that the early mantle may well have been reduced (for 
example, «WM) and that some nodule suites brought to the 
surface in kimberlites and alkali basalts reflect this primitive 
condition: even though the bulk content of volatiles would have 
been much greater, the relative proportions of reduced gas 
species would also have been greater. With time, however, the 
mantle has become progressively oxidized, resulting perhaps 
largely from the preferred loss of hydrogen and carbon. Another 
important factor is that if the upper mantle is geochemicaily 
layered, that is with an upper depleted zone overlying a more 
fertile region* 1 , it is likely that the redox condition of each 
component will differ. Depleted zones result from extensive 
volcanism. so that the regional, stratigraphicai and age of the 
mantle zone being sampled in nodules must be considered. In 
mature mantles underlying cratonic regions of the Earth, a 
mixed sampling of oxidized and reduced components such as 
those present in kimberlites will differ from more fertile mantles 
actively undergoing depletion. Thus, depleted and refractory 
upper portions of mature mantle are likely to be relatively 
reduced. Similarly high level but juvenile and fertile regions of 
the mantle, such as those underlying the ocean floors, are 
relatively enriched in volatile constituents and are, therefore, 
more oxidized. The plausible ranges of f 0 , are WM and FMQ, 
respectively. 

The mantle (and crustal) f 0) model that we propose is illus- 
trated in Fig, 6. Five typical tectonic cross-sections are outlined 
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Ft*. * Selected but typical croas- 
MCI torn through the cruel and 
upper marnle illustrating the distri- 
bution of rock types and profiles 
of oxygen lugacilies The highest 
fot is represented by eondiliont 
lets than magnetite-haematite 
(MH) equilibria. With increasing 
reduction the equilibria are NNO 
iNi-NiO). FMQ (fsyalite- 
magnetiu-quaru), WM (wfistite- 
magnetite), and IW (iron-wustjta). 

The perched and depleted litho- 
spheric layer results from sustained 
volcanisni and varies in thick nets 
as a function of mantle maturity 
and tectonic setting. Nodule* from 
the depleted layer arise from 
era ton, rift, and palaeo-tpreading 
environments. Nodule* from tl vt 
fertile aathenosphere are moat 
common in kimberlite* but are also 
present in smaller concentrations 
in continental rift settings, and arc 
only rarely encountered in pataeo- 
spreading terrains. This distribu- 
tion is related to the depth of 
magma generation which is illus- 
trated. In the text we conclude that the reduced state of the mantle that it inferred by some ' ' * results from selective sampling of depleted 
lithospheres. The lower portion of the upper mantle, in the fertile isthenospheres is oxidized and conditions are at FMQ ranging down to 
WM. to Wl. and to iron at the lower mantle-outer core interface. 



and the oxygen fugacitics we have designated are from data 
presented here and cited earlier I - ,, * ,,J *-”" M - 4 *. For the upper- 
most regions of each cross-section, oxygen fugadties range from 
below magnetite-haematite (MH) to FMQ, and representative 
examples are crstonic crust to ocean floor basalts, respectively. 
Oxygen fugadties decrease within the crust and we propose 
that a redox discontinuity is present at the Moho, transforming 
from FMQ to WM. Lower crustal granuiites are most typically 
fcldspathie, with aimandine-rich garnet and Na-Fe 1 *-bearing 
dinopyroxene; hydrous minerals, such as mica and amphibole 
may also be present. Hence these rocks are more oxidized than 
the underlying and depleted durites and harzburgites. Although 
we suggest that the depleted lithosphere is at WM, it could 
conceivably be lower than WM but not much less, because 
spinels in these rocks are iron(in)-bearing M . The fertile 
aathenosphere, below the depleted lithosphere, signals a return 
to more oxidizing conditions (between FMQ and WM), rep- 
resented by Iherzolites and the genesis of highly volatile kimber- 
lites and carbonatites. With increasing depth the redox state is 
WM, followed by WI and in the lower mantle-outer core 
boundary by iron. This model is applicable to four of the five 
cross-sections (Fig. 6), with the sub-Moho region of the subduc- 
tion profile being more complex, resulting from admixtures of 
various redox components. Voicanism arising from there zones 
are mafic but also add-intermediate in composition and are 
more highly oxidized than basalts at spreading axes or undersat- 
urated volcanic* in continental rift systems. We draw attention 
to two' Anal points: (1) the depleted lithosphere varies enor- 
mously in thickness; and (2) only three tectonic settings give 
rise to mantle nodule suites: cratonic, continental rifts, and 
paiaeo-spreading centres. The common factor is that sustained 
basaltic voicanism has produced a depleted tithospheric mantle 
horizon. Magma genesis in cratonic settings is usually at depths 
>150 km, and hence both fertile Iherzolites and depleted 
dunites and harzburgites are entrained on eruption. On the 
other hand, in alkali-rift environments, magmas are generated 
at much higher levels and the opportunity for sampling fertile 
Iherzolites is diminished. Hence, these volcanics most com- 
monly contain a depleted and refractory nodule suite. Melting 
takes place in the redox zone of WM-FMQ and the resulting 
volcanics are FMQ (Fig. 6), sintiiar to those of kimberlites and 
carbonatites (Fig. 3). Oxygen fugacity assessments, therefore, 
which are based on mantle-derived suites will depend on 


whether the nodules arise from depleted tithospheric mantle 
or from the fertile aslhcnosphere. Herein may lie the source 
of the disparity reported in mantle redox conditions. 

Arculus and Delano 4 * have presented comprehensive and 
convincing arguments for the oxygen fugacity for the eariy 
Earth that includes core formation and the origin of the atmos- 
phere. We are in general agreement with most aspects of their 
model but we differ in interpretation. They conclude that the 
mantle is buffered by IW and that the observed state of effusive 
voicanism at FMQ results from oxidation during upward migra- 
tion of magma. The premise for IW is based on data derived 
exclusively from the electrolytic celt technique, and an explana- 
tion to reconcile these low /o, values is auto-reduction from 
particulate graphite 47 , fluid inclusions comprising reduced car- 
bon spedea (for example, CH< or CO), or some property related 
to thermal relaxation of high-pressure minerals when- heated 
to elevated temperatures at 1 atm. However, the most probable 
reason Is that the data are not an artefact of the method used 
but rather that the peridotites they measured (Fig. 3) arc ail 
highly depleted continental rift nodules. Their model is, there- 
fore, based only on the sandwiched and depleted lithosphere 
layer between the Moho and the fertile asthenosphere. They 
have not considered the latter. 

Iron(lii) and (tl) equilibria in natural silicate liquids 4 *, and 
the pressure dependency of oxygen fugacity in magmas, show 
that fot increases with increasing P (as in Fig. 4) and decreasing 
silica activity 4 *, which implies that the f 0 , of magmas 
(specifically mafic) at depth is greater than that at the surface, 
assuming isochemical ascent. The situation is dearly not this 
simple based on the speciation of iron (Fe l4 /Fe l4 ), and iron- 
oxygen coordination in silicate melts as functions of P, T, fa, 
and bulk composition’^ 1 . In aikaii-rich systems (such as 
NaAlSijO*-NaFe 1 *SijO«, or NajO-Fe^O-SiOj), the ratio of 
Fe’Vtotal iron decreases systematically from l to lOkbar, 
plateaus between 10 and 20 kbar, and increases again from 20 
to 30 kbar. Taken in a broader context, and recognizing that 
the experimental compositions are not directly applicable, it is 
nonetheless noteworthy, that the parabolic behaviour of Fe J * 
as a function of P is consistent with the reduced state at low 
pressures determined by Arculus and Delano 4 *, and supports 
the data we present at higher pressures in kimberlites. The 
interval between 10 and 20 kbar is relatively reduced and 
corresponds approximately to that horizon within the iitho- 
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sphere that w* propoM is similarly depleted in volstilcs (Fig. 
6), and hcnca, is mort reduced. 


Conclusions 

Upper-mantle derived ilmenite nodule tuites from kimberlites 
in West Africa are moderately oxidized and lie within redox 
states consistent with the determinations obtained on other 
nodule tuites in kimberlites by thermodynamic and intrinsic /o, 
methods. These data contrast with lower pressure nodule tuites 
in peridotites from alkali basalts which yield highly reduced 
values. A single all-encompassing redox state for the mantle 
cannot be uniquely inferred. Oxidation and reduction of the 
Earth's upper mantle has evolved with lime, and an expression 
of this evolution is provided by the relative thicknesses of redox 
layers ranging from those in cratonic regimes to those at active 
oceanic spreading centres (Fig. 6). Crustal setting and tectonic 
style, along with nodules and voleanism provide the important 
dues to mantle redox conditions. Mantle source regions for all 
but subduction-related voleanism is WM-FMQ. Subducted ter- 
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rains are more oxidised and lie between NN0 (Ni-NiO) and 
MH. A knowledge of the stratigraphies! position that these 
samples represent is fundamental to resolving the redox state 
of the upper mantle. The issue of mantle fugstitics will not be 
settled until more precise data are available permitting ipedfic 
sample horizons to be definitively located. Cretaceous kimber- 
lites and carbonatites were unquestionably generated in an 
oxidized upper mantle below a depleted lithosphere. However, 
whether extreme depletion of such a zone could have led to 
IW* 4 remains uncertain and problematic. A perched reduced 
zone appears probable in our model but an IW condition for 
the entire upper mantle appears unlikely. 
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Structure of mouse kallikrein gene family 

suggests a role in 

specific processing of biologically active peptides 
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The glandular kallikrein gene family comprises 25-30 highly homologous genes that encode specific proteases involved 
in the processing of biologically active peptides. In the mouse all the members of this family are closely linked on 
chromosome 7. The 9.5-kilobase nucleotide sequence of a mouse genomic clone contains one complete kallikrein gene 
fmGK-U, which is expressed in the male mouse suhmaxiltary gland, and the 3' end of another ( , mGK-2A Differences 
in the coding potential of these genes and the amino acid cequences of other known kaltikreins seem to be functionally 
related to the substrate specificity of the different enzymes. 
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APPENDIX 2 

MAGNETIC SPINELS AND RELATED PHASES FROM 
MANTLE-DERIVED KIMBERLITES: ANALYTICAL DATA 
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Table 2 Eoldu hiaberlitic aplasia 


ID-61-DP-U 


Croundaaaa aplaala 



2.5u - 

2.5u 

- 2.56 

- 2.5)i 

- Edge 

Canter 

- 2.5u 

T10 2 

13.58 

1.73 

3.09 

2.96 

TTu 

5.25 

5.41 

2 tO z 

0.20 

0.18 

0.05 

0.19 

0.04 

0.00 

0.03 

kl l 0 3 

2.86 

0.09 

0.02 

0.01 

0.06 

8.32 

8.40 

Cr z O, 

18.01 

1.07 

0.39 

0.36 

0.29 

45.67 

45.48 

Fa 2 ° 3 

25.84 

65.06 

64.40 

64.58 

66.85 

8.58 

8.13 

FeO 

27,62 

31.17 

27.57 

27.80 

30.66 

19.47 

19.85 

MgO 

10.26 

0.76 

3.00 

2.86 

0.71 

11.59 

11.44 

MnO 

0.74 

0.51 

1.73 

1.72 

0.36 

0.46 

0.39 

NIO 

0.15 

0.12 

0.14 

0.14 

0.13 

0.16 

0.19 

Hb 2 0 5 

Q.0Q 

0.00 

0.00 

0.00 

0.06 

0.00 

.0.00 

Total 

99.26 

100.69 

100.39 

100.62 

100.28 

99.50 

99.32 




CATIONS TO 32 OXYGENS 



T1 

2.833 . 

0.395 

0.695 

0.664 

0.257 

1.048- 

1.082 

Zr 

0.027 

0.026 

0.008 

0.028 

0.006 

0.000 

0.003 

A1 

0.933 

0.032 

0.006 

0.005 

0.021 

2.604 

2.634 

Cr 

3.950 

0.257 

0.093 

0.086 

0.071 

9.587 

9.569 

Fe 3 * 

5.396 

14.869 

14.496 

14.525 

15.359 

1.714 

1.627 

Fe 2+ 

6.410 

7.917 

6.896 

6.950 

7.829 

4.324 

4.418 

Mg 

4.242 

0.344 

1.336 

1.273 

0.324 

4.536 

4.539 

Hr 

0.174 

0.130 

0.437 

0.437 

0.094 

0.104 

0.088 

Hi 

0.035 

0.029 

0.034 

0.033 

0.032 

0.034 

0.040 

Nb 

0.G<10 

0.000 

0.000 

0.000 

0.008 

0.000 

0.000 

Total 

24.002 

23.999 

24.001 

24.001 

24.001 

24.001 

24.000 


Each column represents a single analysis. 


24.001 24.000 






















Table 3 Koidu kiaberlitic grouadmaea ilaanltae 


KD-S1-DP-2B 



Edge 

Getter 

Edge 

Edge 

Center 

Edge 

Center 

Iii) 2 

49.89 

50.59 

50.44 

50.30 

49.86 

50.56 

52.18 

ZrOj 

0.11 

0.11 

0.08 

0.19 

0.15 

0.13 

0.00 

Al 2 ° 3 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.02 

cr 2 0 3 

0.58 

1.09 

0.78 

0.35 

0.36 

0.83 

0.64 

F *2°3 

1.92 

2.09 

3.02 

1.10 

1.77 

2.91 

1.41 

?«0 

40.93 

40.74 

41.42 

41.96 

41.69 

41.13 ' 

41.64 

MgO 

1. 32 

2.54 

1.96 

1.94 

2.04 

2.02 

2.40 

MnO 

1.65 

1.13 

1.06 

1.54 

1.42 

1.09 

1.03 

NiO 

0.23 

0.07 

0.05 

0.10 

0.07 

0.03 

0.07 

Nb^ 

1.07 

0.87 

0.60 

1.60 

1.77 

0.31 

0.11 

Total 

98.20 

99.23 

99.41 

99.08 

99.13 

99.01 

99.50 



CATIONS TO 

3 OXYGENS 




Ti 

0.956 

0.953 

0.952 

0.956 

0,947 

0.957 

0.978 

Zr 

0.001 

0.001 

0.001 

0.002 

0.002 

0.002 

0.000 

A1 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.001 

Cr 

0.012 

0.022 

0.016 

0.007 

0.007 

0.016 

0.013 

f. 3+ 

0.037 

0.039 

0.057 

0.021 

0.034 

0.055 

0.026 

F* 2+ 

0.872 

0.854 

0.870 

0.887 

0.881 

0.866 

0.869 

Mg 

0.069 

0.095 

0.074 

0.073 

0.077 

0.076 

0.089 

ito 

0.036 

0.024 

0.022 

0.033 

0.030 

0.023 

0.022 

Mi 

0.005 

0.001 

0.001 

0.002 

0.001 

0.001 

0.001 

Mb 

0.012 

0.010 

0.007 

0.018 

0.020 

0.003 

0.001 

Total 

2.000 

1.999 

2.00Q 

1.999 

1.999 

1. 999 

2.000 

Mole? 








Hem 

1.86 

1.96 

2.87 

1.05 

1.69 

2.77 

1.31 

Ila 

87.59 

86.05 

8';. 48 

88.39 

87.66 

87.25 

87.51 

Geik 

6.93 

9.57 

7.44 

7.27 

7.66 

7.66 

8.96 

Pyroph. 

3.62 

2.42 

2.21 

3.29 

2.99 

2.32 

2.22 


Each, co loam represents a single analysis. 

liea*hemacice; Ila"ilraenite; Geik-geikielite; Pyroph»pyrophanite. 


Table 3 . Koidu kiaberlitic grcundaass ilaenites 


0-81-DP-22 



BPS 

Edge 

sTTm 

Cenear 

Edga__ 

WB . 

Cancer 

Canter 

Tifl 2 


305 “ 

50.95 

50728 


' '49.94 

Zr0 2 

0.04 

0.03 

0.18 

0.19 

0.11 

0.09 

0.08 

Al 2 ° 3 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

Cr 2°3 

1.00 

0.35 

0.47 

0.51 

1.10 

0.55 

0.88 

P* 2 ° 3 

3.15 

2.55 

1.47 

2.69 

3.21 

2.50 

2.62 

FeO 

41.12 

39.93 

40.31 

40.98 

41.39 

41.14 

40.52 

MgO 

2. 28 

3.25 

3.15 

2.58 

1.95 

2.04 

2.24 

MnO 

1.15 

1.24 

1.14 

1.04 

0.97 

1.16 

1.12 

NiO 

0.08 

0.09 

0.09 

0.07 

0.10 

0.06 

0.07 

Ub 2 0 5 

Total 

0.55 

0.86 

1.53 

0.72 

0.tI3 

0.99 

0.71 

100.52 

99.59 

98.86 

99.73 

99.74 

98.45 

98.18 

Ti 

0.951 

0.958 

CATIONS 

0.953 

TO 3 OXIGENS 
0.955 0.947 

0.952 

0.953 

Zr 

0.001 

0.000 

0.002 

0.002 

0.001 

0.001 

0.001 

Al 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

Cr 

0.020 

0.007 

0.009 

0.010 

0.022 

0.011 

0.018 

f* 3+ 

0.059 

0.048 

O.Qti 

0.051 

0.061 

0.048 

0.050 

f. 2+ 

Q.854 

0.830 

0.846 

0.855 

0,867 

0.873 

0.860 

u g 

0.084 

0.120 

0.118 

0.096 

0.073 

0.077 

0.085 

Mn 

0.024 

0.026 

0.024 

0.022 

0.021 

0.025 

0.024 

Ni 

0.002 

0.002 

0.002 

0.001 

0.002 

0.001 

0.001 

Mb 

0.006 

0.010 

0.017 

0.008 

0.007 

0.011 

0.008 

Total 

2.001 

2.001 

1.999 

2.000 

2.001 

1.999 

2.000 

Mole* 

Has 

2.98 

2.40 

1.40 

2.55 

3.08 

2.40 

2.52 

I In 

86.13 

83.00 

84.43 

85.63 

87.44 

87.39 

86.52 

Geik 

8.47 

12.00 

11.78 

9.61 

7.36 

7.71 

8.55 

Pyroph. 

2.42 

2.60 

2.40 

2.20 

2.12 

2.50 

2.41 


Each, column represents a single analysis. 

Eem-heaatita; Ilm-ilmenite; Geik«*geikielite; Pyroph-pyropnanica, 





Tab 3 Kc .a kiaberlitic groundness oxides 


0-81-DP-2B 


Hasnite# ___ 


1.70 


51.66 
20 
02 
86 


99.97 100.26 100.47 

CATIONS TO 3 OXYGENS 


-In 

49.19 


100.78 


Cmfi p ^S3^t Ctatir Cintir 
i.aF 5.14 4.91 5.02 


69 47.47 49.60 45.29 

42 8.34 8.95 8.96 

19.03 19.21 18.86 

11.94 11.48 11.96 

.59 
.27 
1 

99.08 100.27 100.64 100.01 
CATIONS TO 32 OXYGENS 


Urn 87.02 


Pyropn. 3.12 


Each colunn represents a single analysis except for Ten. ilaenite 
Xen-xenoerystic ilaenite Coverage cf -4 analysis) ; Hsn«heaatite; Ha 
ilaenite; Geifc-geikielite ; ?yToph“pyrophanite. 











Table 3 Kotdu klabarlicic groun idmaaa spinels 


tD-61-DP-2B 

tio 2 

1 

4, *81 

2 

5 .*28 

3 

4? 93 

4 

5 ."ll 

5 

5741 

6 

4.15 

Mt. 

0.09 

ZrO z 

0.10 

0.14 

0.06 

0.13 

0.33 

0.10 

0.09 

ai 2 o 3 

8.82 

8.18 

8.74 

8.74 

5.95 

8.05 

0.00 

Ctj0 3 

48.76 

46.57 

46.36 

45.83 

47.56 

47.18 

0.05 

F *2°3 

8.89 

8.51 

8.88 

8.63 

9.69 

8.68 

67.89 

FcO 

18.58 

19.53 

18.97 

13.85 

13.83 

18.67 

29.71 

MgO 

12.30 

11.82 

11.95 

12.03 

12.17 

12.22 

0.31 

MnO 

0.42 

0.60 

0.53 

0.68 

0,31 

0.51 

0.01 

Nio 

0.24 

0.17 

0.24 

0.23 

0.06 

0.17 

0.67 

il6 2°5 

0.04 

0.06 

0.02 

0.08 

0.13 

0.11 

0.05 

Total 

100.96 

100.86 

100.68 

100.31 

100.94 

100.64 

98.87 




CATIONS TO 32 

OXYGENS 



ti 

0.941 

1.040 

0.969 

1.009 

1.075 

0.976 

0.020 

Zr 

0.013 

0.018 

0.008 

0.016 

0.042 

0.013 

0.013 

A1 

2.708 

2.529 

2.697 

2.705 

1.853 

2.488 

0.000 

Cr 

9.628 

9.655 

9.593 

9.511 

9.939 

9.780 

0.011 

Fe^ 

1.743 

1.680 

1.749 

1.705 

1.927 

1.713 

15.901 

Fe 2 * 

4.048 

4.283 

4.151 

4.138 

4.162 

4.093 

7.734 

Mg 

4.774 

4.619 

4.663 

4.707 

4.793 

4.775 

0.144 

Mn 

0.092 

0.134 

0.117 

0.151 

0.181 

0.113 

0.002 

HI 

0.050 

0.036 

0.051 

0.048 

0.012 

0.035 

0.167 

Hb 

0.005 

0.007 

0.003 

0.010 

0.016 ♦ 

0.014 

0.007 

Total 

24.002 

'24.001 

24.001 

24.000 

24.000 

24.000 

23.999 


Each column rep resent 3 a single analysis, except Mt. 
Mt*magn*tite in serpentinized olivine Caverage of 3 analysis) . 











Table k Koldu fciaberlitic xenocrystlc ilmenites 


D-81-DP-3B 



1 

2 

3, 

4 

CaaeerfS) EdpefS) 

Very 

Bd«e^ 3 ' 

I10 2 

50.38 

51.41 

52.10 

50.30 

49.02 

48. 11 

49.01 

Zz0 2 

0.17 

0.17 

0.14 

0.23 

0.10 

0.01 

0.12 

ai 2 0 3 

0.32 

0.23 

0.17 

0.24 

0.39 

0.31 

0.23 

Cr 2°3 

3.07 

3.90 

3.07 

3.97 

1.91 

3.66 

2.07 

y «2°3 

■ 10.53 

9.37 

9.88 

13.18 

14.74 

16.82 

16.45 

F«0 

18.58 

19.80 

17.67 

16.78 

18.91 

16.62 

16.31 

MgO 

14.77 

14.56 

16.02 

15.62 

14.03 

14.61 

15.30 

MaO 

0.60 

0.65 

0.82 

0.68 

0.61 

0.56 

0.68 

N±0 

0.05 

0.15 

0.05 

0.11 

0.07 

0.16 

0.04 

Kb 2 0 3 

0.16 

0.21 

0.18 

0.07 

0.43 

0.12 

0.17 

Total 

98.63 

100.45 

100.10 

101.18 

100.21 

100.98 

100.38 

T± 

0.872 

0.877 

CATIONS 

0.882 

TO 3 OXTGENS 
0.847 0.843 

0.818 

0.834 

Zr 

0<0Q2 

0.002 

0.Q02 

0.002 

0.001 

0.000 

0.001 

A1 

0.009 

0.0Q6 

0.004 

0.006 

0.010 

0.008 

0.006 

Cr 

0.056 

0.070 

0.055 

0.070 

0.035 

0.066 

0.037 

Fe 3 * 

0.182 

0.160 

0.167 

0.222 

0.254 

0.286 

0.280 

f* 2+ 

0.358 

0.376 

0.333 

0.314 

0.362 

0.314 

0,309 

Hg 

0.507 

0.492 

0.538 

0.522 

0.478 

0.493 

0.516 

Mb 

0.012 

0.012 

0.016 

0.013 

0.012 

0.011 

0.013 

Hi 

0.001 

0.003 

0.001 

0.002 

0.001 

0.003 

0.001 

Nb 

0.002 

0.002 

0.002 

0.001 

0.004 

0.001 

0.002 

Total 

2.001 

2.000 

2.000 

1.999 

2.000 

2.000 

1.999 

MoleS 

Hem 

17.19 

15.38 

15.84 

20.73 

22.97 

25.91 

25.04 

ria 

33.81 

36.15 

31.59 

29.32 

32.73 

28.44 

27.64 

Geik. 

47.88 

47.31 

51.04 

48.74 

43.22 

44.66 

46,15 

Fropfi. 

1.13 

1.15 

1.52 

1.21 

1.08 

1.00 

1.16 


Each, colunn represents a single analysis. 

Kea-fiematite; Hm«ilaenite; Geli-gelltlellte; Pyroph-pyrophanice. 




Tabic 4 • Eotdu kioberlitic spinels. 


KD-81-DP-3B 


Spinal Handing 11a 05 Groundmass Spinal 


Ti0 2 

19.79 

27.23 

18.47 

18.07 

Cancer 

1 U ? 

Edge 
l6 "46 

Zr0 2 

0.05 

0.17 

0.11 

0.03 

0.00 

0,31 

Al 2 ° 3 

3.26 

1.56 

3.06 

3.27 

3.32 

7.13 

Cr 2°3 

3.56 

3.18 

5.97 

3.81 

14.15 

2.14 

P *2°3 

31,78 

16.13 

28.61 

29.05 

27.57 

42.54 

reO 

28.51 

42.31 

30.65 

39.80 

29.43 

29.58 

MgO 

13.60 

8.52 

10.95 

4.78 

9.82 

7,37 

HnO 

1.04 

1.28 . 

0.92 

1.25 

0,91 

0.55 

N10 

0.13 

0.22 

0.10 

0.12 

0.15 

0.20 

NbjOj 

0.00 

• °- Q0 

0.00 

0.00 

0.00 

0.00 

Total 

101.72 

1Q0.60 

98.84 

100.18 

100.13 

100. 78 



CATIONS TO 

32 0XTGENS 



Ti 

3.943 

' 5.687 

3.849 

3.886 

3.061 

2.164 

Zr 

0.007 

0.023 

0.014 

0.005 

0.000 

0.042 

Al 

1.017 

0.509 

1.000 

1.103 

1.078 

2.313 

Cr 

0.746 

0.698 

1.309 

0.861 

3.082 

0.4 66 

Fa 3 * 

6.337 

3.373 

5,966 

6.254 

5.717 

8.808 

f. 2+ 

6.319 

9.831 

7.102 

9.521 

6,784 

6.808 

Mg 

5.372 

3.528 

4.523 

2.039 

4.033 

3.228 

Ma 

0.232 

0.301 

0.215 

0.304 

0.212 

0.128 

Ni 

0.027 

0.050 

0.023 

0.027 

0.032 

0.043 

Mb 

0.000 

0.000 

0.000 

0.000 • 

0.000 

0.000 

Total 

24,000 

24.000 

24.001 

24.000 

23.999 

24.000 


Each coluan represents a single analysis. 





Table 4 Koidukinberlltic spin els. 


O— 81— DP-3B 


Croundaass spliel* 



Cora 

Cora 

Mantle 

Cora 

Mantle 

I10 2 

6.38 

6.11 

18.26 

4.43 

14.69 

Zr0 2 

0.13 

0.15 

0.15 

0.16 

0.11 

ai 2 0 3 

4.23 

5.82 

3.10 

7.62 

3.73 

** 

45.32 

45.90 

28.38 

46.93 

20.58 

? «2°3 

10.62 

10.40 

8.81 

8.92 

22.51 

F«0 

20.67 

20,67 

28.68 

19.61 

27.11 

HgO 

10.85 

11.36 

12,88 

11.00 

11.83 

MnO 

0.74 

0.59 

0.93 

0,39 

0.71 

too 

0.06 

0.12 

0.06 

0.13 

0.15 

Nb,0<- 

0.00 

0.00 

0.04 

0.00 

0.01 

Total 

99.00 

101.12 

101.29 

99.19 

101.43 



CATIONS 

TO 32 QXTGQfS 



Ti 

1.312 ' 

1.219 

3.639 

0.894 

2.954 

Zr 

0.018 

0.020 

0.019 

0,021 

0.015 

Al 

1.363 

1.819 

0.968 

2.411 

1.177 

Cr 

9.795 

9.628 

5.946 

9.958 

4.352 

Fa 3 * 

2.184 

2.076 

1. 757 

1.802 

4.530 

f. 2+ 

4.726 

4.588 

6.357 

4.401 

6.063 

Mg 

4.419 

4.491 

5.089 . 

4.398 

4,714 

Ms 

0.172 

0.134 

0.208 

0.089 

0.161 

Ml 

0.013 

0.026 

0.012 

0.027 

0.033 

Mb 

0.000 

0.000 

0.005 

0.000 

0.001 

Total 

24.0Q2 

24.001 

24.000 

24.001 

24.000 


Each, column represents a single analysis. 



Table 5 b£du xaaoeryecic 1-La tales 


XD-81-DZ8-A2 



1 

2 

3 

4 

5 

6 

Average 

Ti0 2 

48.22 

48740 

48.57 

48.49 

48.86 

39.01 

42792 

Zr0 2 

0.49 

0.60 

0.62 

0.45 

0.58 

0.34 

0.51 

ai 2 o 3 

0.34 

0.02 

0.01 

0,03 

0.01 

0.27 

0.11 

Cr 2°j 

0.41 

0.85 

0.52 

0.54 

0.43 

1.19 

0.66 

? * 2°3 

11.61 

10.55 

12.56 

13.13 

12.59 

29.58 

15.06 

reO 

33.36 

33.70 

27.52 

30.01 

27.47 

21.49 

28.89 

MgO 

5.81 

5.75 

9.21 

7.56 

9.58 

7.77 

7.63 

MaO 

0.20 

0.25 

0.41 

0.33 

0.34 

0.39 

0.32 

»10 

0.00 

0.04 

0.00 

0.05 

- 0.00 

0.00 

0.01 

NbjO, 

0.26 

0.35 

0.30 

0.01 

0.57 

0.44 

0.32 

Total 

100.70 

100.51 

99.72 

100.60 

100.43 

100.48 

100.43 

Ti 

0.876 

0.882 

CATtONS TO 3 OXYGENS 
0.870 0.871 0.86e 

0.706 

0.845 

Zr 

0.006 

0.007 

0.007 

0.005 

0.007 

0.004 

0.006 

A1 

0.010 

0.001 

0.000 

0.001 

0.000 

0.008 

0.003 

Cr 

0.Q08 

0.016 

0.010 

0.010 

0.008 

0.023 

0.012 

f. 3+ 

0.211 

0.193 

0.225 

0.236 

0.224 

0.536 

0.271 

Fe 2+ 

0.674 

0.683 

0.549 

0.600 

0.543 

0.433 

0.579 

Mg 

0.209 

0.208 

0.327 

0.269 

0.337 

0.279 

0.272 

Mn 

0,004 

0.005 

0.008 

0.007 

0.007 

0.008 

0.006 

Ni 

0.000 

0.001 

0.000 

0.001 

0,000 

0,000 

0.000 

Mb 

0.003 

0.004 

0.003 

0.000 

0.006 

0.005 

0.003 

Total 

2.001 

2.000 

1.999 

2.000 

2.000 

2.002 

1.997 

MolaZ 

Rea 

19.22 

17.72 

20.29 

21.22 

20.16 

42.68 

24.02 

ria 

61.38 

62.72 

49.50 

53.96 

43.87 

34.47 

51.33 

Geik. 

19.03 

19.10 

29.49 

24.19 

30.33 

22.21 

24,11 

Pyroph 

0,36 

0.46 

0.72 

0.63 

0.63 

0.64 

0.53 


Each, column represents a single analysis except for the average. 
Kes-hematice*, Ila-iloenice; GeUwgeiiielice; ?7roph*pyTOphanlce. 















Table 5 Koidu kiaberlitic oxides 


i ^fv^TV^e -7 v ’ " T 


E3-81-OZB-A2 

Spinal lamellae Man ding Hamlets on xenocryatic Iloanit7 

Ceone *dj 


tio 2 

8 

2ff?90 

Average 

g.gV" 

1 

49713 

2 

47788 

3 

48720 

49796 

5 

44783 

z«) 2 

0.20 

0.27 

0.26 

0.26 

0.25 

0.22 

0.38 

a 1 2 0 3 

1.13 

0.38 

0.10 

0.00 

0.06 

0.04 

0.13 

Cr 2 0 3 

2.27 

2.18 

2.54 

0.22 

3.11 

1.59 

3.22 

F *2°3 

25.85 

17.96 

12.94 

5.02 

14.66 

12.16 

20.78 

PeO 

45.97 

49. 2Q 

25.14 

27.70 

22.94 

24.36 

20.07 

MgO 

2.81 

3.52 

10.48 

0,31 

10.82 

11.32 

11.13 

MnO 

0.19 

. 0.20 

0.53 

16.75 

. 1.23 

0.54 

6.70 

NiO 

0.20 

0.J.4 

0.01 

0.09 

0.12 

0.02 

0.02 

*2 0 S 

0.03 

0.17 

0.05 

1.96 

0.09 

0.06 

0.11 

Total 

99.55 

99.79 

101.18 

100.19 

101.48 

100.27 

101.37 


CATIONS TO 32 

omras 

CATIONS TO 

3 0X5TGENS 



Ti 

4.634 

5.551 

0.858 

0.913 

0.839 

0.875 

0.781 

Zr 

0.028 

0.039 

0.003 

0.003 

0.003 

0.003 

0.004 

A1 

0.394 

0.302 

0.003 

0.000 

0.002 

0.001 

0.004 

Cr 

0.529 

0.503 

0.047 

0.004 

0.057 

0,029 

0.059 

Fa 3+ 

5.738 

3.949 

0.226 

0.096 

0.255 

0.213 

0.362 

f . 2+ 

11.340 

12.021 

Q.489 

0.588 

0.444 

0.474 

0.389 

Mg 

1.236 

1.531 

0.363 

0.012 

0.373 , 

0.393 

0.385 

tin 

0.047 

0.049 

0.010 

0.360 

0.024 

0.011 

0.014 

Mi 

0.048 

0.034 

0.000 

0.002 

0.002 

0.000 

0.000 

Mb 

0.005 

0.022 

0.000 

0.022 

0.001 

0.001 

0.001 

Total 

23.999 

24.001 

1.999 

2.000 

2.000 

2.000 

1.999 

MoleZ 








Hem 



20.77 

9.09 

23.27 

19,52 

31.48 

Ilm 



44.94 

55.68 

40.51 

43.45 

33.83 

Geik 



33.36 

1.14 

34.03 

36.02 

33.48 

Pyropa 



0.92 

34.09 

2.19 

1. 01 

1.22 


Each column represents a single analysis, except for the average 
Hem-heaatite; Ilm-ilaenite; Geik-geikielita; Pyroph-pyrophanite. 













Table 5 

Koldu fciaberlitic ilaanicas 


XD-61-DZB-A2 

ED-81-DZS-B 





Grouadaass llaanitas 




1 

2 

3 

4 

5 

6 

xio 2 

51727 

49701 

48780 

52748 

50739 

47742 

Zr0 2 

1.03 

0.72 

0.85 

0.34 

0.55 

0.69 

A1,0, 

0.00 

0.00 

0.00 

0.00 

0.01 

0.10 

Cc 2°J 

0.15 

0.14 

0.27 

0.25 

0.04 

0.48 

F *2°3 

7.50 

11.50 

12.94 

9.54 

6.75 

11.03 

FeD 

16.94 

21.15 

22.43 

13.21 

21.36 

24.01 

MgO 

16.77 

12.73 

11.78 

15.00 

10.61 

10.74 

MaO 

2.81 

2.40 

2.21 

4.05 

7.45 

2.07 

HiO 

0.00 

0.00 

0.08 

0.00 

0.02 

0.00 

Nb ? 0, 

2.77 

1.65 

1.29 

1.55 

2.06 

2.05 

Total 

99.24 

99.30 

100.65 

101.42 

99.24 

98.59 

Ti 

0.880 

CATIONS TO 
0.863 0.854 

3 OXYGENS 
0.890 

0.900 

0.853 

Zr 

0.011 

0.008 

0.010 

0.004 

0.006 

0.008 

A1 

0.000 

0.000 

0.000 

0.000 

0.000 

0.003 

Cr 

0.003 

0.003 

0.005 

0.004 

0.001 

0.009 

7. 3+ 

0.129 

0.203 

0.227 

0.162 

0.121 

0,199 

f. 2+ 

0.323 

0.414 

0.437 

0.343 

0.424 

0.481 

Mg 

0.571 

0.444 

0.409 

0.504 

0.376 

0.383 

Mn 

0.054 

0.048 

0.044 

0.077 

0.150 

0.042 

si 

0.000 

0.000 

0.002 

0.000 

0.000 

0.000 

Sb 

0.029 

0.017 

0.014 

0.016 

0.022 

0.022 

Total 

'2.000 

2.000 

* 2.002 

2.000 

2.000 

2.000 

MoleS 

■Tptti 

11.98 

18.30 

2Q.32 

14.92 

11.30 

13.01 

Ha 

29.99 

37.33 

39.12 

31.53 

39.59 

43.53 

Geik 

53,02 

40.04 

36.62 

46.41 

35.11 

34.66 

Pyroph 

5.01 

4.33 

3.94 

7.09 

14.01 

3.80 


Each caiman represents a single analysis. 


Hem-hematite; ria-ilaaaice; Geifc^geikielita; Pyroph-pyrophaaice. 


Table 5 Koldu kiaBarlitic oai du. 
E0-81-DZ2-A2 


Mineral "H" 



8 

9 

10 

11 

12 

13 

14 

TtOj 

17713 

lC68 

11741 

1?T26 

17765 

13.87 

16.97 

ZrOj 

0.2Q 

0.30 

0.18 

0.35 

0.20 

0.26 

0.43 


1.79 

1.64 

2.44 

2.13 

1.60 

1.88 

1.90 

CtjOj 

3.77 

2.27 

13. 12 

1.38 

2.19 

8.29 

1.95 

FeO 

64.85 

66.72 

60.65 

68.76 

65.07 

63.46 

66.02 

MgO 

11.03 

10.27 • 

9.65 

9.93 

10.60 

10.12 

10.17 

MnO 

0.99 

1.08 

0.96 

0.91 

0.98 

0.85 

0.92 

NiO 

0.08 

0..1Q 

0.09 

0.14 

Q.13 

0.12 

0.05 

MbjOy 

0.00 

0.00 

0.00 

0.00 

•0.00 ‘ 

0.00 

0.00 

Tocal 

99.64 

99.06 

98.50 

99.86 

98.42 

98.85 

98.41 




CATIONS TO 

7 OXTGDIS 



Ti 

0.859 

0.852 

0.582 

0.828 

0.899 

0.707 

0.869 

Zr 

Q. 007 

0.010 

0.QQ6 

0.012 

0.007 

0.009 

0.014 

A1 

0.140 ‘ 

0.131 

0.195 

0.170 

0.128 

0.150 

0.152 

Cr 

0.199 

0.122 

0.7Q4 

0..074 

0,117 

Q.444 

0.105 

f* 2+ 

3.605 

3.790 

3.441 

3.894 

3.686 

3.599 

3.760 

Mg 

1.096 

1.040 

0.976 

1.002 

1.070 

1.023 

1.032 

Mn 

0.056 

0.062 

0.055 

0.052 

0.056 

0.049 

0.053 

Ni 

0.004 

0.005 

0.005 

0.008 

0.007 

0.007 

0.003 

Hb 

0.050 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

Total 

5.966 

6.012 

5.964 

6.040 

5.970 

5.988 

5.988 


Each- column represents a single analysis 














Table 5 Xoldu kiabarlitic oxides 


B>-r8l-DP-DZB-A2 


Mineral "R" Spinal Cora* 



• U 

16 

17 

18 

19 

1 

2 

T10 2 

13772 

14770 

itfTT? 

itfTTl 

17723 

5792 

5737 

Zro 2 

Q.15 

0.22 

0.41 

0.69 

0.39 

0.08 

i! 0.22 

*1j0 3 

2.19 

2.23 

2.01 

1.90 

1.98 

3.28 

3.89 

cr 2 0 3 

12.09 

11.29 

2.55 

1.14 

2.21 

47.57 

46.97 

? *2°3 

— 

— 

— 

— 


11.67 

12.34 

FeO 

60.19 

59.48 

67.21 

67.49 

66.44 

21.28 

20.50 

WgO 

11.24 

11.22 

10.05 

9.87 

10.83 

10.53 

10.90 

MaO 

0.99 

0.98 

0.91 

0.89 

0.95 

‘ 0.73 

0.79 

HIO 

0.11 

0.12 

0.14 

0.08 

0.09 

0.11 

0.03 

®2®5 

0.00 

0.00 

0.00 

0.00 

0.00 

0.03 

0.00 

Total 

100.78 

100.04 

100.17 

98.97 

99.12 

101.31 

101. 21 



CATIONS TO 

7 OXTGENS 



CATIONS TO 
32 QXIGENS 

Ti 

0.679 

0.715 

0.851 

0.866 

0.825 

1.200 

1.125 

Zr 

O.QOS 

0.007 

0.013 

0.023 

0.Q13 

0.010 

0.028 

A1 

0.169 

0.172 

0.159 

0.152 

0.157 

1.045 

1.231 

Cr 

0.624 

0.585 

0.135 

0.061 

0.118 

10.157 

9.970 

Fe 3+ 

— 



— 

— 

— 

2.367 

2.494 

f. 2+ 

3.288 

3.262 

3.767 

3.844 

3.759 

4.796 

4.604 

Mg 

1.094 

1.097 

1.004 

1.002 

1.092 

4.231 

4.364 

Mn 

0.055 

0.055 

0.052 

0.051 

0.055 

0.166 

0,179 

Ml 

0.006 

0,007 

0.007 

0.005 

0.005 

0.023 

0.006 

Mb 

0.000 

0.000 

0.000 

0.000 

0.000 

0.003 

0.000 

Total 

5.;ro 

5.900 

5.988 

6.004 

6.024 

'13. 998 

24.001 


Each, colunn represents a single analysis. 


I 


Table. 5 Eoldu klabaglicic «piaala 
. 0*-41 tDZB-A1 


Oroandaaaa tp inala 

Cantar 2 3 4 Edaa Cantar Mid41a 

T10- " 07 " 0751 0.08 0731 11.46 0.07 0.09 


FajOj 69.29 67.81 


local 101.48 101.10 101.68 101.86 101.8Q 101.80 101.12 101.37 

CATIONS TO 32 OXYGENS 


0 


15.815 15.512 15.793 15.706 .10.401 15.308 15.310 15.259 












Table 5 Koidu kiaberlitic spinels 


KD'-31-DZ3-A2 




Spinel Cores (Coat’d) 




3 

4 

5 

6 

7 

8 

9 

T10 2 

6756 

6748 

7.58 

6714 

5759 

7709 

6741 

Zr0 2 

0.35 

0.11 

0.08 

0.12 

0.16 

0.13 

0.13 

ai 2 0 3 

0.41 

3.40 

3.59 

2.47 

3.17 

2.06 

3.06 

Cr 2°3 

2.58 

44.50 

38.82 

46.94 

47.20 

41.32 

46.88 

j Vj 

54.80 

13.68 

14.16 

• 12.01 

12.45 

15.48 

11.10 

FeO 

28.35 

21.23 

23.86 

20.92 

20,66 

21.60 

21.01 

Mg.-O 

4.93 

10.80 

9.10 

10.47 

10.60 

10.24 

10.75 

MnO 

1.33 

0.81 

0.79 

0.88 

0.81 

0.95 

0.77 

NiO 

0.03 

0.04 

0.10 

0.05 

0.00 

0.04 

0.04 

Nb, 0* 

0.07 

0.00 

0.01 

0.00 

0.00 

0.00 

0.00 

Total 

99.41 

101.05 

98.09 

100.00 

100.64 

98.91 

100.17 




CATIONS 

TO 32 OXYGENS 



T£ 

1.460 

1.313 

1.598 

1.265 

1.139 

1.482 

1.310 

Zr 

0.050 

0.014 

0.011 

0.016 

0.021 

0.018 

0.013 

A1 

0.144 

1.080 

1.185 

0.798 

1.014 

0.677 

0.987 

Cr 

0.605 

9.490 

8.607 

10.165 

10.123 

9.084 

10.085 

f« 3+ 

12.203 

2.776 

2.987 

2.475 

2.542 

3.239 

2.272 

f. 2+ 

7.015 

4.789 

5.596 

4.793 

4.686 

5.023 

4.782 

Mg 

2.174 

4.344 

3.803 

4.273 

4.286 

4.245 

4.360 

Mn 

0.333 

0.185 

0.189 

0.205 

0.187 

0.224 

0.177 

Mi 

0.007 

0.009 

0,024 

0.011 

0.001 

0.008 

0.008 

Mb 

0.010 

0.000 

0.001 

0.000 

0.000 

0.000 

0.000 

Total 

24.001 

24.000 

24*001 

24.001 

23.999 

24.000 

23.999 


Each coluan represents a single analysis. 












Table 5 . Kotdu klaherlitic oxides 


O-81-DZB-Al 





Area I 






B 

C 

Z 

A 

E 

F 

D 

a 

X10 2 

17.57 

13731 

2765 

13788 

14773 

6777 

18756 

16752 

Zr0 2 

0.11 

0.08 

0.15 

0.25 

0.15 

0.09 

0.22 

0.27 

^3 

1.95 

2.13 

3.49 

2.17 

2.11 

3.42 

1.78 

1,99 

Ct 2 0 3 

0.95 

11.35 

55.22 

5.15 

9.83 

44.12 

2.48 

2.99 

r *2°3 

— 

— 

8.12 

39.92 

— 

12.61 

— 

— 

FeO 

55.80 

50.55 

19.11 

29.63 

60.57 

21.66 

63.20 

64.57 

MgO 

10.21 

10.21 

9.72 

9.25 

10.58 

10.53 

10.97. 

10.71 

MnO 

1.03 

0.98 

0.67 

0.96 

1.01 

0,76 

1.02 

0.90 

NiO 

0.01 

0.01 

0.06 

0.10 

0.04 

0.03 

0.10 

0.05 

Mb 9 0^ 

0.00 

0.00 

o!oo 

0.00 

0.QQ 

0.00 

0.00 

0.00 

Total 

98.75 

98.62 

99.20 

101.31 

99.02 

99.99 

98.37 

98.00 

Cations 

To 

Oxygens 

7 

7 

32 

32 

7 

32 

7 

7 

T1 

0.901 

0.673 

0.552 

2.890 

0.738 

1.388 

0.935 

0.845 

Zr 

0.004 

0.003 

0.019 

0.034 

0.005 

0.012 

0.007 

0.009 

A1 

0.157 

0.169 

1.134 

0,707 

0.166 

1.099 

0.141 

0.160 

Cr 

0.052 

0.604 

12.033 

1.127 

0.518 

9.512 

0.138 

0.161 

F. 3 * 

f. 2+ 



1.685 

4.406 

8.318 

6.861 





3.787 

3.408 

3.377 

4.939 

3.544 

3.672 

Mg 

1.032 

1.024 

3.995 

3.816 

1.051 

4.280 

1.090 

1.086 

Ma 

0.059 

0.056 

0.157 

0.225 

0.057 

0.175 

0.058 

0.052 

SI 

0.000 

0,000 

0.013 

0.022 

0.002 

0.006 

0.006 

0.003 

Nb 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

Total 

5.992 

5.937 

23.999 

24.000 

5.914 . 

23.999 

5.919 

5.988 


Each column represents a single analysis. 


. Table. 5 Kbidu Klaberlitic spinels 






. KD-81-DZB-B 

KD-81-DZB-A1 



Center 

Mantle 

Center 

Mantlal2 

Center 

Middle 

Mantle 

TiD 2 

5.91 

15.76 

5.77 

50.70 

2.80 

5.66 

14.89 

Zr0 2 

0.13 

0.16 

0.02 

0.30 

0.01 

0.07 

0.24 

A1 2 0 3 

3.03 

1.97 

3.07 

0.00 

3.81 

3.35 

2.30 

Cr 2°3 

47.32 

1.94 

46.80 

0.22 

55.90 

46.29 

1.54 

F. 2 , 3 

11.37 

39.37 

11.35 , 

8.48 

7.45 

12.77 

— 

FeO 

21.12 

30.48 

21.69 

23.27 

18.67 

20.67 

69.92 

MgO 

10.26 

9.50 

9.69 

12.46 

10.17 

10.45 

9.93 

MnO 

G.84 

0.95 

0.70 

2.18 

0.64 

0.79 

0.83 

MiO 

0.14 

0.17 

0.13 

0.00 

0.06 

0.15 

0.10 

NB,0, 

0.00 

0.00 

0.00 

1.79 

0.00 

0.00 

0.00 

Total 

100.12 

100.30 

99.22 

99.40 

99.51 

100.20 

99.75 

Cations 

To 

Oxygens 

32 

32 

32 

3 

32 

32 

7 

Ti 

1.214 

3.307 

1.200 

0.892 

0.575 

1.159 

0.948 

Zr 

0.017 

0.022 

0.003 

0.003 

0.002 

0.009 

0.007 

A1 

0.975 

0.648 

1.002 

0.000 

1.229 

1.075 

0.173 

Cr 

10.225 

0.428 

10.232 

0.004 

12.085 

9.970 

0.078 

Fe 3+ 

2.338 

8.267 

2.362 

0.149 

1.532 

2.618 

— 

Fe 2+ 

4.827 

7.112 

5.016 

0.455 

4.270 

4.710 

3.722 

Mg 

4.179 

3.953 

3.994 

0.434 

4.144 

4.243 

0.942 

Ma 

0.195 

0.224 

0.165 

0.043 

0.149 

0.183 

0.045 

III 

0.030 

0.039 

0.028 

0.000 

0.013 

0.032 

0.005 

Mb 

0.000 

0.000 

0.000 

0.019 

0.000 

0.000 

0.000 

Total 

MoleS 

Sen 

Tim 

Geik 

Fyroph 

24.000 

24.000 

24.002 

1.999 

13,78 
42.09 
40,15 
3.. 98 

23.999 

23.999 

5.920 


Each column represents a single analysis. 

Hem-hematite; Ilm-ilaenite; Geik-geikielite; Pyroph-pyrophanite. 
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ABSTRACT 

MAGSAT (MGST 4/81) anomalies are similar over the West African 
and Amazonian eratons. Differences between total field (aB) and 
N-S horizontal vector (aX) anomalies over both eratons are largely 
those that would be expected for anomalies due to induced magnet- 
ization. At the present low magnetic inclinations, the negative 
AX anomalies which are spatially associated with the cratonic shields 
imply normally-directed positive magnetization-contrast sources . 
Archean nucleii of the shields appear to be more strongly magnet- 
ized than Proterozoic segments, but both stand out from adjacent 
basins and fold belts. Under specified conditions, the uncertainty 
inherent in the poorly defined MAGSAT anomaly zero level nay possibly 
be able to be reduced by co -analysis of anomaly profiles taken 
along opposite sides of the continental rejoin in a pre-drift re- 
construction. Synthesis model anomalies generated from litholog- 
ical and gravity data reproduce reasonably well the cratonic neg- 
ative MAGSAT anomalies. The dominant contribution to MAGSAT AX 
anomalies over the West African and Amazonian eratons results from 
magnetization contrasts within the crust. These are modulated 
by the effects of structural contrasts, intrusions and metamorph- 
ism within the lithosphere, and are localized by the preferred 
distribution of iron deposits into the more ancient crustal blocks. 


lit**" 
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INTRODUCTION 

Since the Introduction of satellite-measured global-scale total field 
magnetic anomaly (&B) maps from POGO (Regan et al . . 1975) and MAGSAT (Langel 
et al .. 1982b) a variety of techniques have been employed In Interpretation 
of the anomalies, including; 1) visual comparison of anomalies to large scale 
tectonic features of the Earth (e.g., Frey , 1982), 2) comparison of anomalies 
to upward-continued aeromagnetlc anomalies (e.g., Lange! et al ., 1980a; Taylor 
et a! ., 1981; von Frese et al . . 1982b; Mon and Son , 1982), 3) comparison of 
magnetization anomalies calculated by equivalent source inversion techniques 
( Mayhew , 1979) to heat flow ( May hew, 1982b), to crustal thickness ( Schnetzler 
and Allenby , 1983), and to the geological characteristics of continental crustal 
blocks (e.g., Mayhew et al. . 1980; Mayhew , 1982a; Mayhew and Gall 1 her, 1982; 
Lange! and Thornlnq, 1982), and 4) comparison of magnetic anomalies to gravity 
anomalies (e.g., von Frese et al .. 1981a; 1981b; 1982a; Taylor et al .. 1981). 
Although Haggerty (1978; 1979) has suggested a possible contrib- 
ution from high Curie temperature metal alloys, most analyses 
have employed the simplyfying assumption, due to Wasilewski et al. 
(1979) » that the mantle is non-magnetic, -and, as summarized by 
Taylor et al. (1983) and by Harrison (1983), have concluded that 
the anomalies may be the combined result of surface geological 
features, lateral contrasts in intracrustal magnetization, and 
variations in thickness of a magnetized layer due either to 
magnetite Curie isotherm (580°c) topography or to crustal 
thickness . 

There are still some outstanding problems, however. Relatively 
high magnetization intensities (1 - 10 A.m."^ - ) have been inferred 
from several lines of evidence (e.g., Regan and Marsh . 1982; 
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Wasllewaki and May hew . 19821 McLeod . 1983 » Schnetzler and Allenby . 
1985). The MAGSAT anomalies are relative to an arbitrary and varying base 
level ( Langel et al . . 1982b; 1982c). Little use has been made to 
date of the MAGSAT vector anomaly maps produced by Langel et al . 
(1982c), except for the work by Galliher and Mayhew (1982). And, 
although on a global scale continental regions tend to appear as 
the most strongly magnetized regions ( Regan et al .. 1975 » Frey . 
1982), most of the quantitative analyses have been made within a 
single continent. 

In this paper we apply simple geological criteria in a qual- 
itative evaluation of the initial (MGST 4/81) MAGSAT maps ( Langel 
et al . . 1982bj 1982c) over West Africa and northern South Americas 
more recent global anomaly maps confirm the reality of most of the 
features of the MGST 4/81 maps ( Cain et ail . . 1934} Langel and Frey . 
1982 ). 

Our strongest interest is in the Man shield of the West African craton 
and in its contrasts and relationships to the surrounding areas as portrayed 
by MAGSAT* The shield contains abundant kimberlites, whose distribution and 
included megacrysts and xenoliths provide information on regional structural 
treads and petrology, geochemistry and mineralogy (e.g., Haggerty , 1932; 

Haggerty , 1933; Haggerty and Tompkins , 1933; Haggerty and Toft , 1934) of 
significance to geophysical models of the local lithosphere. 

Our regional framework for analysis of the MAGSAT aap3 incorporates aspects 
of lithology, gravity anomalies, and the Earth's main magnetic field. 

This is not an evaluation of the MAGSAT aaps solely in relation to continental 
drift. Sather, it is an attempt to address some outstanding problems related to 
the interpretation of satellite magnetic data, and i3 particularly an effort to 
analyze portions of the MAGSAT vector anomaly maps frcm the perspective of a 
synthesis of geological and geophysical data. 
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REGIONAL GEOLOGY 


We examine first the regional geology as ground truth data. 

The West African and Amazonian cratons, which were stabilized before 
1500 m.y. ago ( Clifford. 1970; de Almeida et al . . 1981), are shown In Figure 
1 which outlines the areas enlarged In subsequent Figures where data are 
compiled on the Bullard et al . (1965, fig. 5) 500 fathom (» 1000 m) physio- 
graphic reconstruction of the pre-drift fit of Africa and South America. The 
fit has survived well ( Smith and Hal lam, 1970; Sclater et al ., 1977; Rabi nov/1 12 
and LaBreoue , 1979; Martin et al. , 1981), and for hand compilation and visual 
examination of geological and geophysical features away from the continental 
margins has the advantage that It shows latitude and longitude for both 
continents in a roughly equal-angle, equal-area configuration. 

The generalized regional geoTogy of West Africa and northern South America 
is shown In Figure 2 and locations are referred to Figure 3. Marginal to the 
cratons are fold belts composed of rocks with much younger radiometric age 
dates. In the range 1100 - 350 m.y. ago ( Hurley et al ., 1967). Within both 
cratons the shield areas, uplifted regions of exposed metamorphic and igneous 
rocks, are primarily of two whole-rock age groups - an Archean group (a 2500 
m.y. old) and a Lower Proterozoic group (2200 - 1800 m.y. old) ( Hurley and 
Rand, 1973). Basin or platform deposits partially obscure the craton and fold 
belt structures, but contribute to the lithological details of the reconstruction, 
the principal points of which were enumerated by de Almeida and Black (1972) 
and most recently reviewed by Torquato and Cordani (1981), 
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Shield areas 

Both era tons display a large scale basin and swell structure, with the 
* 

Guiana and Guapore shields north and south respectively of the Amazon basin 
and the Regulbat and Man shields north and south of the Taoudenl basin In West 
Africa (Figures 2 and 3). 

Archean exposures are more vast In West Africa than In South America; both 
are characterized as mlgmatlzed granitic and amphibol e-pyroxene-garnet banded 
gneisses, granites, amphibolites, charnockltes and ferruginous quartzites; 
granite Intrusion and the most recent regional metamorphism, to amphibolite 
or granulite facies, occurred about 2800 - 2500 m.y. ago. The Archean and 
Proterozoic zones in the Man and Guiana shields are separated by major 
mylonltlzed fault zones. Rocks of the Proterozoic age group are also characterized 
as gneisses, mlgmatltes and syn- and post-tectonic granites but also Include 
amphibolites, phyllites, schists, and quartzites which are only occasionally 
ferruginous; these rocks were most recently regionally metamorphosed, up to 
greenschist or amphibolite facies, and intruded by synkinematlc granites, at 
about 2200 - 1800 m.y. ago In a major orogenic event termed the Eburnean in 
southern West Africa and the Trans-Amazonlan In South America (for West Africa : 
Louis , 1970; Dillon and Sougy , 1974; Bessoles , 1977; Black , 1980; Choubert et 
al., 1971; for South America: de Almeida et al ., 1973; Choubert, 1974; de 
Almeida et al 1981 ; Gibbs and Barron, 1983; correlative: Hurley et al ., 

1967; Choubert, 1969; de Almeida and Black , 1972; Gro ss , 1973; Hurley and Rand , 
1973; Torquato and Cordani , 1981). 

A third, relatively minor, uplift in the West African craton exposes a 
small area of Proterozoic rocks to the north of the Tindouf basin (Dillon and 


Soughy , 1974). Uplift in the Ahoggar region east of the craton has exposed 
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2700 and 2000 m.y. old rocks (Hurley et al .. 1967), broadly similar in lithology 
to those in the West African craton ( Black and Girod . 1970)* 

Other African shield regions decrease in age along the coast of 
Africa to 1300 - 800 m.y. old in the Namaqualand shield of the 
Kalahari craton ( Clifford , 1970» Hurley and Rand , 1973)* 

In the Guiana shield of the Amazonian craton the Archean age of the rocks 
south of the Roralma sediments (Figures 2 and 3) Is still open to question 
( Bosma et al ., 1983), and rock and mineral ages younger than 1800 m.y. old are 
common in the west and south of the shield ( Kroonenberg , 1982; Priem et al . , 
1982); most of the Guapore shield Is not older than Trans-Amazonian. Else- 
where In northern South America the S5o Fransisco craton was extensively 
rejuvenated during the Interval 1500 - 1100 m.y. ago and subsequently at about 
500 m.y. ago ( de Almeida et al ., 1973; 1981; Mernick. 1981). Finally, the 
$3o Luis craton in northern Brazil lithologically ( de Almeida et al ., 1981) 
and radiochronologically ( Hurley et al. . 1967) closely resembles the Eburnean 
segment of the Man shield and lies adjacent to it in the Bullard reconstruction 
( Torquato and Cordani, 1981). 

Fold belts 

The fold belts in Figure 2 delineate regions of Precambrian or Paleozoic 
deformation subsequent to the Eburnean - Trans-Amazonian orogeny, and Include 
unmetamorphosed sediments and wide ranges of metamorphic and igneous rock types; 
Pan-African - Brasil iano cycle deformation dated as afcout 450 +_ 250 m.y. old 
( Hurley , 1973) is generally restricted to the metasedimentary margins of the 
fold belts adjacent to craton regions, whereas the gnelssic and qranitic 
infrastructures of the fold belts frequently display Lower Proterozoic and 
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and Archean ages ( de Almeida et al. , 1973; Grant, 1973; Machens, 1973; Bessoles, 
1977; 1980; de Almeida et al. . 1981; Choubert and Faure-Muret, 1971). Suture belts 
In Figure 4, from Frey (1979) and Burke et al. (1977) are shown without respect 
to age, and do not exactly correspond to the fold belts In Figure 2. 

The tectonic origin of these fold belts Is not, necessarily, an open and 
shut case, but the conflicting hypotheses of an enslallc in-sltu origin ( Hurley. 
1973; Kroner, 1977; McEl hinny and McWilliams. 1977; McWilliams and Kroner. 

1981) and a Wilson-cycle tectonics origin ( Burke et al ., 1977; Black et al ., 

1979) may be generally reconciled by assuming the opening and closing of 
marginal, narrow seas. Exposed ophiolitic assemblages are relatively rare, 
but see, for example, Leblanc (1976), de Almeida et al . (1981), Fountain and 
Salisbury (1981) and Berbert et al. (1981). In either case, the fold belts 
are in*. ~->'tant in the regional tectonic control of the locus of Mesozoic southern 
South Atlantic rifting, generally parallel to the fold belt structural trends 
in both continents ( de Almeida et al . , 1973; Nairn and Stehli , 1973; Torouato 
and Cordani , 1981). 

Basina and continental rifting 

The major basins in Figure 2 are characterized by depth to 
basement isopachs ( Bronner et al ., 1980; de Almeida et el ., 1931) indicating 
the areas of greatest subsidence. The Intracratonic basin is older in West 
Africa than in northern South America: deposition within the Taoudeni basin 
was strongly developed between 1100 - 650 m.y. ago, primarily toward a trough 
in the north (see Figure 2) but 4 - 8 km of sediment are also associated with 
the Gourma aulacogen in the south (Dillon and Sougy, 1974; Black et al . , 1979; 
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Bronner et al .. 1980; Lesouer and Houssine-Pouchklne . 1980), whereas the 
sedimentary record In the Amazon basin begins In the Silurian with subsequent 
deposition being directed toward three subsidiary basins separated by the 
Iqultos and Purus arches at about 70°W and 60°W respectively ( Bigarella. 1973; 
de Almeida et al. . 1981). A much earlier cycle of sedimentary deposition Is 
also recorded, In the Roralma Formation on the Guiana shield and the Tarkwaian 
deposits on the eastern Man shield, prior to 1700 m.y. ago ( Hurley and Rand, 

1973). The major part of the northern South American sedimentary record and 
the graben surrounding the Atlantic Ocean (see Figure 4) developed In response 
to stress associated with Mesozoic continental rifting ( Burke and Whiteman . 

1973; Nairn and Stehli . 1973; Burke , 1976; Burke et al . , 1977;). 

Volcanic features 

Injections of basaltic dikes and sills and eruptions of flood basalts were 
also Intimately associated with the separation of the continents (e.g., Bigarella , 
1973; Basaltic Volcanlsm Study Project ., 1981, pp. 33 - 35), as was the slightly 
later kimberlite volcanism evident in the =120 m.y. old dike-swarms and pipes 
of Liberia and Sierra Leone ( Haggerty . 1982); kimberlites, from Dawson (1980) 
and Bardet (1973), are shown in Figure 4 as a guide to potential lower crustal 
and upper mantle sampling sites in relation to large scale magnetic anomalies. 
Carbonatitic intrusions post-date the kimberlites at about 60 m.y. old ( Marsh , 
1973; Moore , 1976; Herz , 1977), and volcanic sites active within the last 1 m.y. 
(Trombka and Lowman, 1979) are included in Figure 4 as a primitive guide to 


high heat flow regions. 
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Nature and distribution of Iron ore deposits 

In the analysis of any magnetic anomaly maps It should be natural first to 
ask whether there are any known highly magnetic bodies in the region in question. 
Iron ore deposits In West Africa and northern South America shown In Figure 4 
are the total resources (l.e., both known reserves and potential ores* United 
Nations (1970)) as of 1970, from Marelle and Abdullah (1970) and Alvarado (1970) 
slightly updated for South America ( Dorr* 1973a; 1973b; Tolbert et al .* 1973). 

Within the major era tons the largest resources are magnetite- and hematite- 
rich oxide and oxide-silicate facies bedded iron-formations of the Lake Superior 
type ( Gross * 1965* 1970). Many of these deposits are referred to as itabirltes, 
i.e., laminated metamorphlc rocks dominantly or entirely composed of granular 
quartz and iron oxides ( Parr, 1952), corresponding in general terms to the 
ferruginous quartzites which are most prominently developed in the 2700 m.y. 
old and, to a lesser extent, the 2000 m.y. old shield rocks ( Marelle and Abdullah, 
1970; Alvaradg 1970; Goodwin, 1973; Gruss, 1973). They were probably originally 
chemical sediments deposited without a strong volcanogenlc association in shallow 
sea or continental shelf type of environments ( Gross, 1970; 1973), and have been 
subsequently metamorphicallv deformed and super- and hypogene enriched into 
substantial Iron ore deposits. 

The massive variety of iron deposit ( Gross , 1970) is well developed in the 
Atlas Mountains ( Marelle ?nd Abdullah, 1970) and In Peru and Chile ( Alvarado, 
1970), originating either from iron carbonate sediments, or magmatic activity, 
for example the El Laco magnetite lava flow in Chile ( Park, 1961 ). Residual 
deposits, the third and last major variety of iron deposit ( Gross, 1965 ; 1970), 
are common as lateritic crusts on basic intrusive rocks throughout the areas of 
wet tropical climate ( Alvarado, 1970; Marelle and Abdullah, 1970). 


GRAVITY ANOMALIES 


Relating magnetic anomalies to gravity anomalies is a 
common approach to analysing magnetic anomaly data. At a 
cratonic or continental scale the available gravity data sets 
are maps of free air and geoid values and Bouguer anomaly maps. 

The most recent global free air gravity anomaly maps (not 
shown) which have been published in a graphic or map form are 
GEM 8 ( Wagner st al .. 197?) and GEM 10 ( Marsh , 1979)- Values were 
derived from ground-based and satellite-measured data, and the 
high degree and order (n > 12) maps show features of half- 
wavelength *10° with amplitudes of a few mgals. Above the West African 
and Amazonian cratons both maps are dominated by negative anomalies 
of peak amplitude -20 mgal, which overlie the centrally located 
intra-cratonic basins. Positive anomalies to the north and south 
slightly overlap the basin regions and are less closely spatially 
related to the shield regions than the positives are to the basins. 
However, the Guiana shield and particularly the Man shield areas 
both correspond to relatively weak (S+12 mgal) but persistent 
positive anomalies . 

In the NASA GSFC global gravity geoid shown by Taylor fit al . 
(1983), which is based on a combination of GEM 10B- and 1° :c 1° 
surface gravity data, the Amazonian craton is swamped by a 
gradient field from the Andes to the Atlantic coast, whereas in 
West Africa the general area of the Man shield is a prominent 
positive geoid feature but other features are poorly resolved. 


Bouguer gravity anomalies (Figure 5) for West Africa ( Slettene etal. . 
1973) and South America ( Wilcox. 1977) are partially from ground-based observa- 
tions* but are supplemented by statistically* topographically and geophysically 
predicted values In the unsurveyed areas ( Brevllle et al . . 1 973 ; Brown et al . . 
1980} a«e also Orlin , 1966). Only selected contours are shown in Figure 5; 
contours were selected by visual examination and were chosen, with constant 
contour intervals as far as possible, so as to emphasize features with wave- 
lengths >200 km, the expected spatial resolution of MAGSAT (see following 
.section). West African craton values range between +20 and -40 mgal, comp- 
ared to Amazonian craton values ranging between about +40 and -80 mgal. 

Some of the Bouguer anomalies are a strong function of topo- 
graphy, in cart because of the prediction methods. An example is 
the strongly negative Andean mountains. East of the Andes are 
several positive anomalies, overlying basins and plains s the 
Amazon basin, for example, is characterized by a series of short 
wavelength mostly positive anomalies . The S-40 mgal anomalies 
over the Amazonian craton shield areas have trends which align 
with superimposed sedimentary deposits ( Wilcox , 1977), such as 
over the Roraima sediments of the Guiana shield, and over the 
volcano -sedimentary deposits overlying the central region of the 
Guapore shield . 

Within. the West African craton are two prominent relatively positive (> -20 
mgal) Bouguer anomalies over the Taoudeni basin* both of which are related to 
the Gourma aulacogen ( Lesquer and Moussine-Pouchkine , 1980). The western 
boundary of the craton is marked by a shift to positive values associated with 
the Mauritanide fold belt ( Chiron. 1974), whereas the craton is bounded in the 


east by a series of short wavelength positive anomalies (Louis, 1970) roughly 
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apparent in Figure 5 as the trend of the -40 mgal contour. The western Archean 
zones of the shield regions have generally negative (s -40 mgal) Bouguer values 
and appear to be undercompensated (Louis, 1970; Behrendt and Hotorson , 1974), 
but the coastal gradient Is clearly compromising the Man shield Bouguer anomaly. 
The eastern segment of the Man shield shows no pronounced features In Figure 5, 
but Is, In detail, a region of short wavelength highs and lows corresponding 
to small grabens and uplifts (Hastings , 1983), whereas the central region of 
the Regulbat shield Is relatively positive and the eastern section Is swamped 
by a large negative anomaly extending over both shield and basin regions. 


MAGSAT ANOMALIES 


The assembly and the expected spatial resolution of. the MAGSAT 
anomaly maps are critical to their geological interpretation. 
Scalar, or total field, anomalies (aB) and north-south horizontal vector 
Mf&malles (aX) from MAGSAT (4/81) maps by Lange! et al . (1982b; 1982c) are 
transposed onto the Bullard reconstruction in Figures 6 and 7; magnetic inclinations 
(1980) are Interpolated from Langel et al . (1982b) and Dobrin (1976, fig. 14.8). 
The anomaly maps were derived from 550 - 350 km altitude orbital profiles by 
subtracting from screened magnetically "quiet" data : 1) spherical harmonic 
model main field values, 2) potential function estimates of external fields, . 
and 3) linear ramps of a few tens of degrees in latitude; statistically out- 
lying values were then discarded and the residual profiles were compiled and 
contoured ( Langel et al .. 1982b; 1982c). 

The resultant "crustal" anomaly maps show features with half-wavelengths 
on the order of 1° -10° and of magnitude sJl5|nT, but considering 1) that the 
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maps are compromised by * poorly defined zero level, were 
contoured on the basis of 2° x 2° block averages at unequalized 
altitudes, and, at equatorial latitudes, include aB and aX noise 
levels of about 1 nT (Langel et al., I982b| 1982c), 2) that there 
are unresolved problems in the spherical harmonic separation of core, 
crustal and ionospheric sources ( Langel and Estes , 1982; Carle and Harrison , 

1982; Klumpar and Greer , 1982; All dredge , 1983a), 3) that the north-south 
orbit and along- track data processing emphasize east-west striking anomalies 
( Langel et ai . , 1982b), particularly in £ B at low latitudes, and 4) considering 
the results of spectral analyses In the spatial (Sailor et al. , 1982; All dredge , 
1983b) and time ( All dredge , 1983c) domains, it is clear that only persistent 
anomalies exceeding about 2 nT In magnitude may be plausibly related to 
lithospheric rocks, and then only to geological units exceeding about 200 km 
diameter. 

Anomalies over the era tons 

Over both the West African and Amazonian cratons, which are presently 
situated within 40° magnetic inclination of the geomagnetic equator (0° inclination) 
are negative MAGSAT AB anomalies where the main field is close to horizontal 
over the.southern shield exposures, positive anomalies over the central parts 
of the cratons at about 15° north Inclination, and negative anomalies at or 
just above 30° inclination (Figures 2 and 6). These same general features are 
also evident In the P0G0 data, becoming more clearly resolved as the data 
continued to be processed from Regan et al . (1975) through Frey et al . (1979) 
and Langel et al . (1982b). 

According to general principles, and as described by Langel et al . (1982c), 
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induced a 6 and A X anomalies generated over bodies of high relative magnetic 
susceptibility at or close to the geomagnetic equator should be negative and 
very similar in form and location, as are the anomalies over the north and 
south of both cratons (Figures 6 and 7). The amplitudes of the A 8 and a X 
anomalies are closely equivalent at the geomagnetic equator, and di- 
verge with increasing magnetic inclination. If induced, the negative 
AX anomalies should have negative aZ (vertical component) and positive aZ 
anomalies to the north and south respectively, as in fact, they do (compare 
Figure 7 to the A Z map t-y Langel et al . (1982c)). Relationships with a Y 
(horizontal east-west) are equally clear in theory, but are less obvious in 
practice, possibly because of the poorer overall quality ( Langel et al. . 1982c) 
of the A Y map. 

But although these cratonlc anomalies are persistent and behave as if they 
are induced, they are not well Isolated from adjacent features, and where the 
main field Is very close to horizontal an induced negative A B or a X anomaly 
is centered over its source, but has smaller positive cusps to the north and. 
south. This can lead to a superimposition of anomalies, as described in 
Central Africa by Hastings ( 1982) , so that, for example, the positive AX 
anomaly over the Taoudenl basin may in part be due to the strongly magnetic 
regions to the south and north. 

Induced magnetic anomalies are, furthermore, a function of magnetic 
inclination, complicating any spatial relationships to other data. However, 
where the main field is close to .horizontal, as in these cratonic regions, AX 
is less sensitive to inclination than is a B; reduced-to-the-pole (RTP) 
anomalies, although the standard techniques are non-ideal at low inclinations 
(Langel and Frey , 1983), have maximum and positive values directly over their 
Ignoring possible overlap, at low magnetic inclinations 


sources . 
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the most significant difference between the RTP AB map of South 
America by Hinge et_al . (1982) , the unreduced AB map (Figure 6), 
and the AX map (Figure 7)» should be in the geographic latitude 
location of a given anomaly peak as a function of magnetic inclin- 
ation. 

This Is tested In Figure 8 where the arcuate line of profile Is selected 
In order to Intersect as nearly as possible the peak values of the profiled 
anomalies In all three maps, and to be roughly perpendicular to magnetic Inclina- 
tion contours. Relative to the RTP anomalies the majority of the profiled 
anomaly peak shifts can reasonably be attributed to the effects of magnetic 
Inclination: the peak displacements are greater for a B than AX, are symmetric 
about the geomagnetic equator, and Increase In magnitude poleward. This relation- 
ship Is also evident over the rest of South America and the surrounding oceanic 

areas, is further support for an induced origin of the anomalies, 
and is subsequently employed in geographically relating the Figure 7 AX 
anomalies to geology and gravity. 

Anomalies related to the shield areas 

In both cratons it is the shield areas which spatially corr- 
espond to negative AX anomalies indicating regions of positive relative 
magnetization. Over the Amazonian craton the Guiana shield negative 
(<-2 nT), which shifts progressively slightly equatorward from 
aB to AX to RTF, corresponds to the veneer of Roraima sediments, 
to a Bouguer low and topographic high, to the Arche an iron-rich 
exposures north of the Roraima, and it closes at the edge of the 
continental shelf. The Guapore shield negative anomaly, at about 
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0° inclination, corresponds to (Phanerozoic) cover, low Bouguer 
values, exposed Proterozoic shield, and is interrupted along 50° 
west longitude where the craton is bounded by a marginal fold belt. 

For West Africa, Hastings (1982) visually associated the 
Archean segments of both West African shields with stronger 
negative AB values than occur over the younger shield segments, 
and Figures 2 and 5 suggest a similar relationship to Bouguer 
anomalies. In respect of the South American inclination relat- 
ionships (Figure 8), the aX anomaly (Figure 7) in the north of 
the West African craton at about 30° - 35° north inclination is 
a saddle with two < -4 nT peaks, seemingly corresponding to the 
Archean and Proterozoic structure of the Reguibat shield. In 
contrast, the Man shield AX anomaly peaks (i - 4 nT) at the geo- 
magnetic equator over an area including Archean shield and 
fold belt features, and decays eastward over the Eburnean Prot- 
erozoic shield segment falling to the 'zero level' by the 
craton boundary. Toward the west, both the Man and Reguibat 
shield anomalies have closures reflecting the shape of the cont- 
inental shelf. 

In geographic comparison to the distribution of iron 
deposits shown in Figure 4, the negative aX anomalies over the 
West African craton have peaks which rather closely correspond 
to the bedded iron-formation deposits of the Archean shield 
segments, particularly bearing in mind the expected latitudinal 
peak shift of the Reguibat shield anomalies. The Guiana 
shield negative AX anomaly may be similarly related to bedded 
iron formation, but the Guapore shield is not known, to date, 
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to have such substantial deposits. Moot generally, however, 
it is simply the shield areas of both cratons which are assoc- 
iated with negative AX MAGSAT anomalies. 

Given that the West African and Amazonian cratonic shield 
areas are all four characterized by similar negative MAGSAT AX 
signatures, how general a result is this ? The Sao Luis craton 
(Figure 1) is radiometrically, lithologically and structurally 
closely related to these two cratons, particularly to the Eburnean 
segment of the the Man shield. Prom the West African - Amazonian 
craton shield-MAGSAT relationship the expected signature of the 
Sao Luis craton would be a weak negative MAGSAT AX anomaly, and 
where the craton meets the coastline ' (Figure 2) is a small (*200 
km) relatively (<0 nT) negative AX peak at 2.5°S 43°W (Figure 7) • 

Even if poorly resolved at the limit of resolution of 
MAGSAT, the anomaly appears to be well related to the SSo 
Luis craton, thus lending support to the combined West 
African and Amazonian craton interpretive model of 
MAGSAT anomalies. However, the model cannot be immediately and 
succesfully applied to the adjacent Sao Fransisco or Zaire cratons, 
and it is not suggested as a global model but rather as particular 
to the presently widely separated West African, Sao Luis and 
Amazonian cratons . 

Anomalies along the continental rejoin 

Our primary interest in the shields is presented here in the framework of 
a pre-drift reconstruction, and in placing POGO and MAGSAT maps of Africa and 


South America adjacent to each other Frev et al . (1983) and Gaideano (1981; 1983) 
have previously noted a general correspondence across the continental rejoin of 
A B anomalies, mostly the positive ones. The rather low degree of correspondence 
In detail, however, can be judged from Figure 6, and It Is appropriate, there- 
fore, briefly to examine anomalies along the continental rejoin. 

Profiles of MAGSAT A X values for the African and South American sides of 
the continental rejoin are presented In Figure 9, which was constructed from 
the AX map with Interpolated values at other than the even Integer contour 
values. Between African latitude 30°S and the point where the longitude line 
10°W intersects the African coast, the line of the profile bisects the overlap 
or underlap (mismatch) of the 500 fathom (* 1000 m) isobath; west of 10°W for 
approximately 500 km the 110 m.y. old oceanic magnetic lineaments, from Sclater 
et al . (1981), are the lines of the profiles, which are thus located, for the 
most part, a few tens of kilometers seaward of the 500 fathom iso- 
bath. Vectors at about 300 km station intervals indicate the (1980) 
apparent magnetic inclinations in the curvilinear planes of the 
profiles and are proportional in length to the (1980) MAGSAT main 
field model values of Langel et al . ( 1980b } 1980c). 

These profiles should not be overinterpreted, either in relation to each 
other or as to spatial coincidences with known structures, because 1) the zero 
nT base level varies in an unknown manner over the global MAGSAT maps ( Langel 
et al ., 1982b; 1982c), 2) the profiles are derived and interpolated from already 
smoothed 2° averages, and 3) are subject to different and rather wide ranges 
in main field inclination and intensity. 


Nevertheless, the profiles give an immediate impression 
similar to that of slightly out-of-phase sine waves relative to 


different base levels . Some earlier studies have suggested 
that certain POGO and MAGSAT anomalies in this area arise from 
structures pre-dating Mesozoic continental drift ( Gald^ano . 1981 > 
Hastings . 1982 » Hinze et al . . 1982 » Frey et al . , 1983 » Galdeano . 
1983)* The combined West African - Amazonian craton interpretive 
model certainly supports this idea in detail. 

Very detailed intercontinental comparisons are compromised , however, 

by the poorly defined nature of the MAGSAT base level, or MAGSAT anomaly 
aero level. It seems possible that, in the future, the 2 ero level ambiguity 
might be able to be reduced by forcing profiles like those in Figure 9 
to agree with each other. Assuming similar characteristics of source volumes 

and magnetizations, geographically opposed FTP anomalies could be co-nor realized 
in amplitude across a continental rejoin, allowing one of the apparent zero 
levels to float. Anomalies from two continents possibly could then be analysed 
with respect to a similar base level. The wide range of magnetic inclinations 
precludes such a use of the Figure 9 profiles. Even with co-normalized HTP 
profiles the absolute base level still would be unknown, but the general method 
may be a useful approach to the zero level problem. 


DISCUSSION 

The previous sections have examined rather separately the 
geology and the gravity and MAGSAT anomalies over West Africa 


and northern South America. • The following discussion attempts to 
synthesize these observations and to pinpoint possible origins 
for the MAGSAT anomalies over the cratons, particularly with 
respect to the Man shield of the West African craton. 


Surficial distribution of iron 

A striking but surficial MAGSAT-geological correlation 
emerging from the West African - Amazonian craton interpretive 
model is the association of iron ore deposits with regions which 
are implied from the MAGSAT anomalies to be strongly magnetic. 

The relationship applies only to certain of the Precambrian 
bedded iron- formations within the two cratons, however, and is 
not suggested as a global or general relationship. Exposed iron- 
formations, including non-economic iron-enriched sequences, do not 
have the volume to be solely responsible for a MAGSAT anomaly, and 
the Liberian (West Africa) itabirites, for example, are drastically 
attenuated by minimal filtering applied to conventional aero- 
magnetic data ( Behrendt and Wotorson , 197 * 0 * 

But in relation to the ultimate objectives of the MAGSAT project ( Lange 1 
et al ., 1982a) the correlation may be (Significant for resource exploration . 
because, although there certainly are exceptions, even a cursory examination 
shows that the correlation occaisionally is apparent outside of the West African 
and Amazonian cratons. 

The area of the Bangui anomaly, for example ( < -12 nT at 5°N 
18°W in Figure 7), includes only minor amounts of iron- formation, 
but for an induced crustal origin the anomaly requires a body 


about 3 x 10^ kro^ of susceptibility contrast at least 0.125 S.I.* 
( Regan and Harsh , 1982) , which approaches that of an iron-formation 
susceptibility contrast to surrounding non-magnetic rocks. A second 
example is the Kursk positive anomaly (MAGSAT AE >+15 nT at 50°N 
18°E); at the local magnetic inclination this intense anomaly 
implies a large positive magnetization-contrast source. The feature 
is dominant in early PCGO maps (e.g. , Regan et al . , 1975» fig* 1) and has 
been frequently but only briefly mentioned in the satellite magnetic literature 
(e.g,, Coles et al . , 1982$ Frey , 1932$ Frey et al . , 1983$ Meyer, et 
al . . 1983$ Olivier et al .. 1983)* Corresponding to the anomaly is 
an extremely iron-rich region, containing in 1970* almost 12 % of 
the world's known iron-ore resources ( Sokolov , 1970$ U nited Nations 
1970, p>5), the protores of which are bedded iron-formations which 
were regionally metamorphosed at about 2700 and 2000 m.y. ago 
( Alexandrov . 1973; Flaks enko et al . , 1973)* 

Despite these coincidences, the distribution of iron-rich 
regions in West Africa and northern South America can be only 
indirectly related to MAGSAT anomalies. It seems probable that 
both phenomena are present-day and surface manifestations of deep- 
seated and long-lived processes which have affected the geochemical 
differentiation of iron, and which may have produced other large 
scale inhomogeneities in the crust. 

Long wavelength aeromagnetic surveys 

We turn now briefly to examine regional aeromagnetic trends 
for information which could contribute to an understanding of 
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MAGSAT anomalies over the West African and Amazonian shield regions. 
Aeromagnetic coverage within- the cratons is too sparse to allow 
compilation and upward continuation for comparison with MAGSAT anomalies, but 
analyses over the Liberian portion of the Man shield ( Behrendt and Wotorson , 

1974} and In Minas Gerais, Brazil ( Gasparini et al .. 1979; 1981), tend to 
confirm the suggestion by Strangway and Vogt (1970) that the longer wavelength 
U 30 km) aeromagnetic anomalies over Precambrian shields in both continents 
are dominated by regional features which are related to the basement structure. 
Behrendt and Wotorson (1974), for example were just able to distinguish Proterozoic 
and Archean shield segments in their long wavelength aeromagnetic data. Green 
(1976) observed similar relationships in Project MAGNET data over central and 
southern Africa, and It Is quite feasible that certain MAGSAT anomalies repre- 
sent the integration and amalgamation of the same deep-seated basement contrasts 
in magnetic properties as are responsible for aeromagnetic anomalies. 

This suggestion Is supported by Wasilewski and Mayhew (1982) who summarized 

lowf.i crustal intensities from a number of long wavelength aeromagnetic studies 

In different regions of the world as most frequently being in the range 1-10 
-1 

A.m. . The aeromagnetic study in the Superior province of the Canadian shield 
by Hall (1974) exemplifies this type of observation, although this particular 
area corresponds in the MAGSAT maps to a weak anomaly south of a 
strong positive anomaly which is associated with exposed high meta- 
morphic grade rocks ( Coles et al . , 1982). Similarly, the lower 
crustal or upper mantle a 5 A.m. sources below the Ukraine shield 
( Krutikhovskava and Pashkevich , 1979s Krutikhovskaya et al ., 1979) 
are swamped by the Kursk anomaly. 

The implication is clearly that lower-crustal high magnetic 
intensity sources suggested by long wavelength aeromagnetic anom- 
alies may be, but are not necessarily always, a cause of satellite- 
measured magnetic anomalies • 

.. - ** h 
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Geological and geophysical structural contrasts 

In all probability it is contrasts in both chemistry and 
structure that are responsible for MAGSAT anomalies. What is 
known about the deep crustal structure of the two cratons, and 
what are the consequences for MAGSAT anomalies ? 

Portions of the Man shield have surface heat flow and heat 
production values ( Sass and Behrendt , 1980) typical of shield values 
throughout the world (c.f., Rao and Jessop , 1975), and there are magnetotelluric 
conductivity contrasts ( Ritz. 1982; 1983) and teleseismic velocity contrasts 
( Briden et al 1981; 1982) compatible with a shield crustal thickness of 
about 40 km. Otherwise, published heat flow or crustal thickness data for the 
West African and Amazonian cratons are very scarce, and it is not possible 
at the present time thoroughly to evaluate the MAGSAT anomalies in 
relation to structurally-based contrasts in magnetization due 
either to Moho or Curie isotherm topography. 

However, considering the fundamental differences normally tc 
be expected between oceanic and continental crust , some of the 
negative AX anomalies, for example the -4 nT anomaly over the Man 
shield, may partly be due to the continental-oceanic contrast in 
crustal thickness . But the Man shield shows contrasts to the 
marginal fold belt in conductivity and seismic velocity structure 
to depths of 200 km ( Ritz , 1982} 1983; Briden et al . . 1981; 1982), 
so that part of the MAGSAT anomaly probably also arises from 
structural contrasts within the continental lithosphere . If the 
deep crustal structure only were known in more detail there would 
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doubtless be other structural factors able to be invoked in explain- 
ing the MAGSAT anomalies . 

Man shield gravity model of MAGSAT anomalies 

The range of observations above widens rather than constrains 
the possible origins of West African and Amazonian craton MAGSAT 
anomalies, so that, although coincident gravity and magnetic 
anomalies need not necessarily be caused by the same source rocks, 
it is qualitatively useful further to examine gravity-magnetic 
relationships over the shield regions, and to ask whether the 
potential field anomalies jointly can be reconciled by a simple 
model . 

As a basis for this there is a persistent and consistent 
mathematically inverse relationship between satellite-derived 
free air and AE magnetic anomalies over the United States, docum- 
ented by von Frese et al . (1982a). Their models are for high 
magnetic inclinations, so that, because of the inclination depend- 
ence of magnetic anomalies, a mathematically direct relationship 
between free air and magnetic anomalies is implied for low magnetic 
inclinations . Neither of the Amazonian shields have definitive 
upward -continued ground-based free air anomalies ( Hinze et al ■ , 
1982); in the GEM 8 and GEM 10 anomaly maps the Man 
shield corresponds to a positive free air and geoid anomaly, but 
this is opposite in sign to the a 3 and AX MAGSAT anomalies. The 
von Frese et al . (1982a) models clearly are not immediately applic- 
able to the West African and Amazonian MAGSAT anomalies. 
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One way in which these observations may be partially recon- 
ciled depends on whether the resolution of density contrasts 
at depth is related to the different scales of the 
gravity surveys. Bronner et al . (1980) interpreted the positive 
Isostatic anomalies of Louis (1970) as due to the combined effects of locally 
thinned upper crust and the relatively excess mass of Iron-formations . The 
positive satellite-derived gravity anomaly which we have associated with the 
Man shield Is from the high degree (> 12) maps, which are believed to represent 
diapiric sources between about 50 - 350 km depth (Kahn, 1971? 

Sclatef et al ., 1975* Marsh and Marsh , 1976). 

The Bouguer anomalies In Figure 5 may originate almost anywhere, from the 
topographic surface on down, however the maximum uncertainty of even 
the predicted values in Africa is + 20 mgal ( Slettene et al. , 1973), insufficient 
to change the sign of the Man shield anomaly, so that if any part is due to a 
feature other than topography the anomaly implies a density deficiency. In 
relation to the Brown and Girdler (1980) interpretation of African gravity, 
particularly the < -100 mgal Ahoggar uplift anomaly, sources for the « -40 mgal 
Man shield anomaly are most reasonably at lower crustal or upper mantle depths. 
The combination of isostatic, free air and Bouguer anomalies suggests a 
crustal thickness of 37 + 7 km according to relationships out- 
lined by Woolard and Strange (1962), in general agreement with the = 40 km thick- 
ness implied by magneto tell uric and teleseismic experiments of Briden et al . 

(1981; 1982) and Ritz (1982; 1983). 

The simplest way in which the geology, gravity and magnetics 
can all be simultaneously reconciled is the following. 
Regional uplift, characteristic also of several other, topographic 
swells in the African continent (LeJBas, 1971; Hi Ison, 1972; Burke and Wilson , 






1972) i and an upwelling mantle source are reflected in the positive 
geoid and free air anomalies j by the Woolard and Strange (1962) 
relationship the Bouguer anomaly reflects the density contrast at 
the seismic Moho $ and the positive isostatic anomaly reflects a 
thinned and iron-rich upper crust ( Bronner et al . , 1980). Assum- 
ing after Waailewski et al . (1979) that the mantle is non-magnetic, 
the seismic Moho is equivalent to the lower limit of magnetization 
in the low heat flow shield area. The negative MAGSAT anomaly 
thus represents the integration of 1) magnetization contrasts within 
upper crustal rocks, 2) magnetization contrasts in the lower crust 
(which was elevated closer to the topographic surface and into 
lower temperature regimes by uplift and consequent erosion), and 
3) the magnetization contrast at the crust-mantle boundary. 

Role of Intrusions and seroentinization in the Man shield model 

If the Moho is not a simple interface of constant 
and all-encompassing character then there may be other possible 
causes for the gravity and magnetic anomalies over the Man shield 
and elsewhere. For example, the lower crustal or upper mantle 
depth (=40 km) density deficiency implied by the Man shield 
Bouguer anomaly, if real, may be further reconciled if alteration 
of the lithosphere due to intrusion and uplift is widespread. 

The Man shield is a locus for Mesozoic and perhaps even 
earlier kimberlite voleanism (Figure 4j Haggerty . 1982) and 
basaltic voleanism also associated with continental rifting 
( Behrendt and Wotorson . 1974). A myriad of very specific 
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tectonic events are generally invoked to assure upper mantle 
magmas an unimpeded path to the surface (e.g., Dawson , 1980), but 
in all probability many of these magmas may never reach the 
surface. The magnetic consequences of solidification at depth are 
twofold i the first is the crystallization of magnetite solid 
solutions in high modal proportions ( Haggerty , 1975) • and the 
second is metasomatism brought about by the release of dissolved 
volatiles, and their infiltration, laterally and vertically. 

Oxides are sensitive to environmental temperature and oxygen 
fugacity, and the upper mantle is considered to be sxrat- 
igraphically layered in oxidation-reduction potential (the 
lower crust and lithosphere are reduced whereas the upper crust and astheno- 
sphere are relatively more oxidized (Haggerty and Tompkins . 1983)). Even though 
the process of serpentinization is one that involves extensive hydration, low 
oxidation states commonly exist ( Moody . 1976), giving rise to magnetite and 
metal alloys (e.g., Haggerty , 1978 and references therein). Whether or not the 
alloys in serpentinltes could contribute significantly to their magnetic proper- 
ties ( Haggerty , 1979; Wasilewski et a! ., 1979), the serpentinization of peridotites 
is frequently a process involving considerable volume expansion ( Moody. 1976), 
and an increase in magnetic susceptibility due to the formation of magnetite, 
concommitant with a decrease in density ( Hatherton, 1967; Saad , 1969; Henkel . 

1976; Lienert and Wasilewski, 1979; Blakely and Page , 1980). 

These factors could be related to the density and magnetic 
contrasts associated with the Man shield gravity and MAGSAT anom- 
alies, but their significance has not been demonstrated on a shield- 
wide scale and is therefore relegated to being just another possible 
contribution to the overall lithologic contrast of the Man shield 
to its adjacent areas. 
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Significance of regional metamorphism 

Perhaps the longest standing contrast related to the West African and 
Amazonian shields stems from their regional metamorphlsms > 1500 m.y. ago. In 
multiply metamorphosed Precambrlan shields, the magnetic oxide minerals of 
gnelsslc and granitic rocks have typically been subject to several episodes 
of reequilibration, and the long range magnetic characteristics of such meta- 
morphic terrains are complex. 

However, the crux of the magnetic anomaly - magnetic 
mineral Interrelationship Is that high T1 contents in mineral oxides decreases 
Curie points, magnetic susceptibility and coercivity, with high magnetic 
Intensities becoming most readily possible upon oxidation, as reviewed by 
Hargraves and Bannerjee (1973). During regional metamorphism titanomagnetite 
in low grade (chlorite zone) rocks may be oxidized toward hematite, ana Fe partition 
into phylloslllcates at intermediate grades ( Rumble , 1976). Titanium from 
surviving titanomagnetite Is concentrated into ''exsol ved" and discrete ilmenite 
and partitioned Into silicate phases at higher grades ( Lidiak . 1974; Rumbl e , 

1976; Hall et al . , 1979). Very fine grained and magnetically stable iron 
oxides may "exsolve" from Iron silicates ( Haggerty , 1979), thus increasing 
magnetization intensity In high grade metamorphlc rocks. 

Magnetizations of granullte facies basic rocks are among the highest values 
measured from a variety of tectonic settings (Wasilewski and Hayhew , 1982), 
and Coles et al. (1982) note stronger MAGSAT and aeromagnetic 
anomalies over high metamorphic grade regions of the Canadian 
shield than exist over adjacent lower grade regions. The strongest 
negative MAC-SAT anomalies over the West African and Amazonian cratons 
appear to correlate with the more ancient and higher metamorphic grade 
Arche an cores of the shields . These complexes are dominated in the 
uuper crust by granitic gneisses, granulites, migmatites and amph- 
ibolites, and their magnetic contrasts are not due to an exotic mineral 
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ogy or to an intrinsic abundance of magnetic minerals, except locally 
where iron-formations are predominant . 

In lower crustal metamorphic granulites from the Man shield, 
however, native iron has been identified, formed by reduction of 
ilnenite and as a decomposition product ■ of almandine-pyrope garnet 
( Haggerty , 1985; Haggerty and Toft . 1984). As with metal alloys due 
to serpientinization, the magnetic significance of native iron in 
lower crustal or upper mantle rocks is possibly only local, but is 
certainly another factor contributing to magnetization contrasts assoc- 
iated with different degrees of metamorphism. 


Numerical simulation of MAGSAT anomalies 


We have suggested above that certain MAGSAT anomalies over West Africa 
and northern South America have sources which are reflected in the exposed 
geology and in gravity anomalies, but such a qualitative visually-based 
interpretation does not fully explain how the magnetic anomalies originate. 

There are fundamental problems in seeking a more quantitative interpret- 
ation, however. Magnetic anomaly sources derived from inverse modelling are 
subject to a lack of unique 'solution, and, as reviewed in the introduction, 
these techniques sometimes give rise to solutions which are themselves anomalous 
in the magnetic intensities required. Forward modelling methods on the other 
hand, require at least some geologically-based geometrical constraints, but, 
els noted above in discussion, geological information on subsurface structure 
is almost completely lacking for the 'West AfricEin and Amazonian cratons. 
Furthermore, forward models require constraints on magnetization and magnet- 
ization contrasts, but the identification cf native iron (Fe°) in lower 
crustal rocks from the Man shield ( IIaggert„ , 1983; Hapjgerty and Toft , 19*34) 
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contraindicates the usual modelling assumption that crustal magnetization 
depends Significantly only on magnetite. 

In view of these problems associated with traditional approaches to 
defining magnetic anomaly sources it is not possible at the present time to 
quantitatively and completely evaluate the sources of MAGSAT anomalies over 
the West African and Amazonian crat-Jiia. However, if there do exist geolog- 
ically meaningful correlations between the MAGSAT anomalies and lithologic, 
metallogenic and gravity features then it should be possible to simulate the 
MAGSAT anomalies based on the hypothesized relationships. 

To attempt this we made the following assumptions and rationalizations: 

1) MAGSAT anomalies represent the integration of magnetization contrasts 
within the upper crust, the lower crust and the upper mantle; 

2) magnetizations within these three gross structural levels contribute 
equally to the anomalous magnetic field at several hundred kilometers altitude; 

3) upper crustal magnetization contrasts are due to lithological contrasts 
exposed at the Earth's surface; 

*0 gravity anomalies are lithologically related to coincident magnetic anomalies 

5) the most significant lower crustal magnetization contrasts are due to 
differences on crustal thickness and are reflected in the Bouguer anomaly field; 

6) other lower crustal magnetization contrasts and contrasts within the upper 
mantle are reflected in free air anomalies; and 

7) determining actual numerical values of magnetization contrasts is less 
significant than the identification of relatively stronger or weaker potential 
magnetic sources. 


following these assumptions, the six areally dominant types of exposed 
geological units of the cratons and their margins (Figures 2 and k) were rank 
ordered in terms of more or’ less strongly magnetic according to observations 
adduced in the earlier stages of discussion and they were assigned small 
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dimensionless ranked numerical values? an allowance for rifts in tectonic 
settings other than a basin was inserted, and the continental shelf was 
assigned a small negative value. The total ranges of Eouguer and GEM 3 
free air gravity anomalies within the cratons were then divided into six 
arbitrary but equal increment ranges and these were also assigned similar small 
dimensionless values (Table 1). For the geological units, it should be 
stressed that the assigned dimensionless numbers are not intended to be absolute 
values of magnetization or even magnetization contrasts, but simply represent 
a rationalized and empirically derived subjective rank ordering of apparent 
relative magnetic intensity. As regards the values assigned to the gravity 
ranges, these numbers are not intended to represent density contrasts directly- 
equated to magnetization contrasts and are not meant to have a physical meaning 
other than to serve as a representation of the assumption that coincident 
gravity and magnetic anomalies are lithologically related. 

To complete the simulation, the cratons and surrounding areas v/ere 
gridded into 2° blocks (^200 x 200 km, corresponding to the spatial reso- 
lution of KAGSAT) and each 2° x 2° 'station' was assigned from Table 1 three 
values according to the dominant type of geology and the gravity anomalies within 
the block; the values at each station were then sunned and the dimensionless 
results v/ere contoured (Figure 10). 

The simulated anomaly nap bears a first order resemblance to the MAGS AT 
anomaly maps (Figures 6 and 7) in the general distribution of the most pronounced 
negative and positive peaks. Actual numerical values of the simulated anom- 
alies are much less significant than the fact that the dominant simulated 
negative anomalies occur over the cratonic shield regions. But the positive 
; IAGSAT anomalies are net so well reproduced by the simulation; fer example , 


in West Africa the Senegal Window and che Gourma aulacogen (Figures 2-4) appear 
as weak negative and positive simulated anomalies, respectively, in place of 
a single positive MAGSAT AB and AX anomaly over the Taoudeni basin, however, 
if the MAGSAT anomaly is caused primarily by an overlap of the positive cusps 
of negative anomalies over shields to the north and south, then the MAGSAT 
positive should not be expected to be reproduced by a simple simulation whfch 
takes no account of the effects of magnetic inclination. 

To account for all features of the MAGSAT anomaly maps would require not 
only a full knowledge of the subsurface geology, and incorporation of the effects 
of magnetic inclination, but also would, require a more precise definition 
throughout of crustal-anomaly versus core magnetic fields. Although not 
quantitative, the simulation has the merit that it is simple and it results 
in ’anomalies' which are reasonably close in character to the MAGSAT anomalies. 
■Ulus the simulation lends strong support to our qualitative interpretation of 
MAGSAT anomalies over the West African and Amazonian cratcns as being due to 
a variety of superposed lithological and geophysical contrasts. 


CONCLUSIONS 


The basic geological relationships between and within the West African 
and Amazonian cratons have long been established and provide a coherent frame- 
work in which to evaluate recently available MAGSAT data. The two cratons 
have a similar geological structure, and the immediate conclusion is that 
the MGST (4/°1) anomaly maps show similar features ever both cratcr.s. 

The behaviour of the MAGSAT anomalies over these cratons is largely consistent 
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with the expected behaviour of* Induced magnetization. A qualitative 
combined analysis of £B and AX maps Indicates that the dominant negative 
anomalies, which must be related at the present low geomagnetic Inclinations to 
sources of positive magnetization contrast and normally-directed 
magnetization, are spatially associated with the Precarabrian 
shield regions of both the West African and the Amazonian cratons. 

Further comparison of A B, AX and RIP MAGSAT anomalies confirms the 
expectation that a X Is less sensitive to inclination effects than is aB. 

The aX anomalies are more closely spatially related to the known upper 
crustal geology In that they more readily resolve the distinction between 
Arehean and Proterozoic shield segments, particularly in West Africa, such 
that the older segments appear to be more strcngly magnetized. A correlative 
example in support of the craton model is the Sao Luis craton, a small wedge 
of Proterozoic shield in northern South America with close affinities to the 

African craton: it is apparent In the AX map as a minor relatively 

negative anomaly. 

In relation to the pre-continental drift reconstruction employed 
to present the various data, profiles of iX constructed along the continental 

rejoin illustrate some of the problems involved in intercontinental 
comparisons of MAGSAT data, and suggest that, under circumstances where certain 
physical parameters of the anomaly sources are known, forcing similar profiles 

to agree v/ith each other may help to establish a more realistic 
base level for the maps. 

Concentrating on the Man shield of the West African craton, the super- 
imposition of isostatic, Bouguer, free air and geoid anomalies 
indicates density contrasts within the upper crust, the lower 
crust and the upper mantle. These can be simply reconciled 
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with the positive magnetization-contrast source implied by the 
negative MAGSAT AX and AB anomalies if the lower crust is at least 
partially, and perhaps substantially, responsible both for the 
Bouguer anomalies and for the MAGSAT anomalies. 

A simulation of MAGSAT anomalies based on the Man shield relationships 
and taking geological and gravity data as input generates 'anomalies' which 
are reasonably close in form to the MAGSAT anomalies over the We st African and 
Amazonian cratons, particularly for the large negative anomalies indicative of 
regions of positive magnetization contrast. The simulation tends to confirm 
that the overriding causes of MAGSAT anomalies in these areas are crustal rocks, 
but the possible contribution of upper mantle sources in providing a background 
magnetization is also accounted as significant. The spatial coincidence of 
several negative MAGSAT anomalies with iron-rich regions of the upper crust 
may be more than just coincidence; in the simulation the iron-formations serve 
to localize the 'anomalies', but in reality both the magnetic anomalies and 
the iron ore deposits are most plausibly surficial manifestations of deep- 
seated geochemical differentiation and regional iron-enrichment. 

This qualitative analysis of MAGSAT maps in relation to a synthesis of 
the geology and geophysics of the West African and Amazonian cratons has shown 
that these tv/o large and now widely separated regions of the Earth are intimately 
portrayed in the MAGSAT vector anomaly maps. That the vector maps are amenable 
to such an interpretation confirms their overall accuracy, consistency and validity. 
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FIGURE CAPTIONS 

Figure 1. The location map identifies the West African (1) 
and the Amazonian (2) cratons, and outlines the areas enlarged 
in subsequent Figures. The other cratons are i 3 = Sao Luis, 4 = 
Sao Fransisco, 5 = Zaire. 

Figure 2. Generalized regional geology (based on de Almeida 
and Black (1972) and Torquato and Cordani (1981) and other ref- 
erences in the text) assembled on the Bullard et al ■ (1965) re- 
construction. The 500 fathom isobath is omitted in areas of 
overlap. Two generations of shield rocks are identified, along 
with fold belts and basins . 


Figure 3« The location map identifies the Man and Reguibat shields of 
the West African craton and the Guyana and' Guapore shields of the Amazonian 
craton, along with other locations referred to in the text (M.F.B. = Mauritanide 
fold belt; G.A. = Gourma aulacogen; S.W. = Senegal window; ROH.. = Boraima 
sediments) and the two major basins, the Taoudeni and the Amazon. 


Figure. 4. The tectono-metallogenic map indicates the major 
graben and sutures ( Burke and Whiteman . 1973 > 3urke . 1976; Burke 
et al . , 1977; Frey . 1979), the Gourma aulacogen ( Lesquer and 
Moussine-Poucnkine , 1980), thrust faults (Bowman, 1979), recently 
active (<1 m.y.) volcanic structures ( Trombka and Lov/man , 1979), 
kimberlites as both dikes and pipes (Bardet, 1973? Dawson . 1980), 
and iron ore resources known in 1970 ( Marelle and Abdullah , 1970; 
Alvorado . 1970). 
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Figure 5. 3ouguer gravity anomalies for West Africa ( Slettene 
et al . . 1973) and South America ( Wilcox . 1976), from ground based 
observations supplemented by predicted values . Shield areas of 
the two major cratons show predominantly strong negative anomalies. 

Figure 6. Scalar, or total field, (aB) MAGSAT (MGST 4/81) 
anomalies at about 400 km altitude ( Langel et al . , 1982b), and 
magnetic inclination contours. At inclinations <35° dominantly 
negative anomalies are induced by positive magnetization contrast 
sources. Shield areas of both West African and Amazonian cratons 
underlie predominantly negative anomalies . The Bangui 'anomaly ( < 
-12 nT) is in Africa at about 15° south inclination. 


Figure 7* North-south horizontal vector (aX) MAGSA.T (MGST 
4/81) anomalies ( Langel et al . , 1982c), which are similar in form 
to the AB anomalies but more closely spatially related to the out- 
cropping shield geology. 

Figure 8 . MAGSAT anomaly profiles across northern South 
America, taken roughly perpendicular to inclination contours and 
through equivalent peaks of the aB and aX maps (Figures 6 and 7) 
and the redueed-to-pole (RTP) A3 map by Hinze et al . (1982). Rel- 
ative to the RTP profile , AX is less sensitive to magnetic inclin- 
ation than is aB. (reversed scale on RT? profile) 


Figure ' 9. MAGSAT AX profiles of A.frica and South America 
taken along the Bullard et a l . (1965) continental rejoin. Arrows 


at each station are proportional in length to the (1980) main 
field strength, and indicate apparent inclinations in the profiled 
planes . 

Figure 10 . Simulated MAGSAT anomalies. Contour values are 
dimensionless . 




TABLE CAPTION 


Table 1 . Input values for numerical simulation of MAGSAT anomalies. 
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TABLE 1 


GEOLOGICAL INPUT 
Bedded iron formation, 


in Archean shield -3 

Archean shield -2 

Proterozoic shield -1 

Continental shelf -0.5 

Fold belt 0 

Sedimentary basin +1 

Rift or aulacogen, 

in basin. +2 

Rift or aulacogen, 
in. other setting add +4 


BOUGUER INPUT 

< -60 mgal .... -3 , 

-60 to -40 2 

-40 to -20 1 

-20 to 0 0 i 

0 to + 20 ... . + 1 
>+20 +2 

FREE AIR INPUT 

< +8 mgal -3 

+8 to +2 -2 

+2 to -4 -1 

-4 to -10 0 

-10 to -16. . .+1 
>-16 +2 
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APPENDIX 4 

NATIVE IRON IN THE CONTINENTAL LOWER CRUST 
PETROLOGICAL AND GEOPHYSICAL IMPLICATIONS 


Abstract. Lower crustal granulite xenoliths recovered from a kimberlite pipe 
in v/est Africa contain native iron (Fe°) as a decomposition product of both 
garnet and ilmenite. Magnetic measurements show that less than 0.1 volume 
per cent of iron metal is present. Geothermometry and oxygen geobarometry 
indicates that the oxide and metal phases equilibrated between iron-wustite 
and magnetite-wustite buffers which may represent the oxidation state of the 
continental lower crust, and the depleted lithospheric upper mantle. Ferro- 
magnetic Fe° could be stable to a depth of = 95 km, and should be considered 
in the interpretation of long wavelength static magnetic anomalies. 


Knowledge of the mineralogy of the deep continental crust is extremely 
limited (1), and the redox state of the lower crust has never been fully 
addressed. Although the Earth's core is probably dominated by metallic iron, 
surface terrestrial conditions are generally oxidized and naturally occurring 
iron metal, native iron, is rarely formed except in coal beds, low temperature 
(<500°C) serpentinites, lavas which have or may have incorporated carbonaceous 
sediments (2) and, at one locality, in a quartz garnet glaucophane lawsonite 
assemblage (3). 

We report the recognition of native iron in lower crustal granulites 
from the Man Shield in Liberia, west Africa. The rocks (2-23 cm diameter) 
were recovered from a diamond-bearing kimberlite pipe (10°4l'W 7 G 33'N) that 
erupted in the Mid-Cretaceous, 90-120 my ago, following the breakup of 
Gondwanaland (4). Granulites and other xenoliths were incorporated from the 
volcanic conduit walls into the kimberlite on its upward passage from the 
mantle. Although high temperatures were attained, rapid adiabatic cooling 
has prevented thermal metamorphism of the xenoliths or chemical interaction 
with the kimberlite. 

The granulites, composed of 1-5 mm size garnet, clinopyroxene and plagioclase 
with accessory ilmenite and rutile, are similar to granulites described from 


elsewhere in Africa (5) and other parts of the world ( 6 ). Garnet-pyroxene 
geothermoaetry (7) yields equilibration temperatures of 700-780 °C assuming a 
10 kb lower crustal pressure, commensurate v/ith temperatures and pressures of 
granulite facies metamorphism. Partial melting in some granulites has 
produced second generation pyroxene with kyanite and scapolite, and these 
rocks record temperatures well above 800°C. Butile is overgrown by 
Fe 2+ -freudenbergite (Ka^FeTiyO^g) which is sodium-metasomatic in origin (8) 
from reaction with a fluid phase also enriched in Ca, S and CO 2 (and possibly 
CO and CH^) that crystallized interstitial scapolite (9) and pervasive sulfides 
(troilite, pentlandite, chalcopyrite) ; rare graphite formed along cracks and 
mineral grain boundaries. 

Garnets are Fe^ + -rich with 51-56 mole fs almandine (Fe~Al-Si_0-_), 27-32 % 

3 2 3 12 

pyrope (Mg^AlgSi^O^) and 14-17 % grossular (Ca^Al^ijO,^) . Incipient to 
advanced decomposition has produced iron metal particles (<2 yin to = 10 ym), 
parti all y or completely enveloped in magnetite. Metal-oxide assemblages are 
concentrated along crystal edges and in irregular fracture patterns throughout 
garnets (Fig. la). Metal compositions, determined by electron microbeam 
analyses, show 98-99 wt % Fe, with minor Mi (=0.15 wt %) , Co (=0.35 wt %) , 

Cu (“0.30 wt J$), and Zn (=0.2 wt ?o). The matrix to native iron is a spongy, 
fine grained, white and translucent phase identified in some cases as 
scapolite, but in other samples it remains unidentified. A highly al u m inous 
( 18-21 wt % Al^Oj) tschermakitic amphibole, of assumed high pressure origin 
(10), is also present. Experiments (11) indicate that almandine garnet is 
unstable at oxygen fugacities ( f 0^ ’ s ) above the FKQ (fayalite = magnetite + 
quartz) buffer at 2 kb. and at temperatures above approximately 850°C. Higher 
pressures increase the thermal stability, however there are no experimental 
data on garnets of the present compositions. 

Native iron also resulted from the decomposition of ilmenite (Fig. 1b). 
The assemblage is ilmenite (FeTi^O^) + ulvospinel (Fe^TiO^) ± iron (Fe°), with 
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troilite (FeS) an accessory phase. The iron has similar grain sizes to, but 

lower Ni, Co, Cu and Zn (total of <0.5 wt ?o) than, particles in garnet. Very 

low values of ferric iron in both the ulvospinels (0.2-0. 5 wt % Fe-0_) and 

2 3 

ilmenites (0.1-2. 4 wt e ,o) are consistent with experimental data (12) for Fe-Ti 

oxides in equilibrium with metallic iron, suggesting that ilmenite decomposition 

took place at =1000°C at low £0^ (<IW, iron-wustite) and that the ulvospinel 

may be cation-deficient and Ti-rich. Applying the solution parameter model for 

metal-free coexisting Fe-Ti oxides (13) temperatures of 996° to 830°C at f0 2 's 
•*12 2 6 2 

between 10” * to 10** * atms are obtained. These data lie just above the 

MW (magnetite-wustite) buffer at 1 bar and lie between MV/ and IV/ if the buffers 

are corrected to 10 kb total pressure. 

From iron metal textural relations in garnet and in ilmenite, and based 

on the differences in temperature obtained from garnet-clinoxpyroxene and 

ilmenite-ulvospinel mineral thermometers, we suggest that native iron formed 

in the solid state at higher temperatures than those implied from Fe-Mg and 

Ca exchange in the silicates. If pressures greater than the assumed 10 kb 

pertained, then higher temperatures would result from the calculation. We 

conclude that decomposition was initiated below IW and equilibrated to IW-MW. 

Native iron was formed at high temperatures ( c 800-1000°C) and decomposition 

was induced at high pressures (=10 kb), based on the assemblage scapolite, 

tschermakitic amphibole and Fe 2+ -freudenbergite. 

With respect to broad scale geophysical surveys of the Earth's crust, 

selected physical properties of the larger specimens of granulites are in 

Table 1. Apparent specific gravities range from 2.75 to 3.25 but are mostly 

between 2.95 and 3.25. Equating specific gravity to density, by Birch's 

relationship for granulites (14) as modified by temperature (15) compressions! 

* 

seismic wave velocities (V ) of 6.3 to 7.2 km. sec - are estimated 

Curie balance thermomagnetic analyses of several specimens confirmed the 
presence of iron metal and magnetite. Each analysis ( = 10 mg chips from samples 
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3 i 4, 9, 25) showed reversible magnetic transitions at 580° and 770 °C (±25°C), 
the Curie temperatures of pure magnetite and pure iron, respectively. 

Magnetic susceptibility (X y ) of the granulites at room temperature and 

-3 

pressure is 0.5-3 x 10 (S.I.), natural remanent magnetization (NRM) 

. . -3 -1 ti 

intensities range from about 10 to 10 A.m , and the Konigsberger ratio is 

usually >1 (Table 1). Figure 2 suggests a proportional relationship between 

the restricted values of susceptibility and the wide range of NEM intensity. 

Long wavelength (50-500 km) aeromagnetic and FOGO and MAGSA? satellite- 
derived magnetic anomalies frequently appear to require lower crustal magnetic 
sources with magnetization intensity contrasts in excess of 1 A.m (17) ♦ a 
range of induced magnetization rarely attained on a large areal scale by 
dominant crustal materials exposed at the Earth's surface (18). Magnetite- 
bearing granulites are considered as the most probable source rocks (19), 
particularly if thermally-promoted viscous magnetic remanence (20) or the 
Hopkinson effect of thermally-enhanced (ex2) susceptibility (21) of magnetite 
are considered, but metal alloys in the upper mantle have also been suggested 
as possible mineral magnetic sources ( 22) • 

Native iron observed in the Liberian lower crustal granulites (Fig. 1) 
is in the multidomain (MD) size range (25), allowing an estimate, based on 
magnetic susceptibility (24), of <0.1 vol % Fe°, assuming that only MD Fe° 
contributes to susceptibility. Even with this same assumption, in accounting 
for the very large Hopkinson effect observed with pure iron (25) the granulites' 
average induced magnetization cannot be enhanced to 1 A.m” or greater until 
depths equivalent to temperatures of 580° to 620°C. Contrasts in total 
magnetization intensity between granulites, however, range from 0.01 to 10 A.m 
at 25°C, and the highest intensities are due largely to NRM (Table 1). 

Pure KD iron does not maintain a directionally stable magnetic remanence 
because the coercive force is “80 A.m“ , close to the magnitude of the Earth's 
main field. Furthermore, the critical temperature of stable remanence 
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(blocking temperature) in heating is usually a spectrum extending below the 
Curie point in MD materials. Also, magnetic minerals can acquire a viscous 
remanent magnetization (VPM) parallel to an external field, particularly if 
exposed for long periods of time, even at temperatures below the blocking 
temperature. Thus, MD Fe° remanent magnetization in the lower continental 
crust could he largely aligned with the Earth's field, and would then be 
indistinguishable from Induced magnetization as far as high altitude magnetic 
surveys are concerned. 

We conclude that these native iron-bearing granulites are examples of 
the types of rocks which may be responsible for a significant part of the 
magnetisation measured by POGO and MAGSAT satellite magnetometers, and that 
the lower continental crust may be characterized by oxidation states equivalent 
to IV7-MW, Unless there is an intermediate relatively oxidized horizon, these 
oxidation states probably persist into the upper mantle and down through the 
depleted lithosphere ( 26 ), suggesting that native iron may also exist below 
the Mono. If so, then temperature and the oFe-yFe phase boundary rather than 
the ct-ct'Fe Curie transition are the significant controls on the lower 
depth limit for iron ferromagnetism in the lithosphere. The ct-y boundary ( 27 ) 
intersects the typical continental shield geothermal gradient at 75 to 80 km, 
but the deepest limit may be defined by the very low gradient Sierra Nevada 
geotherm which intersects the a-yFe phase boundary at SO to 95 km. 

The granulites described here contain metallic iron, which, along with 
the sulfides and graphite but mitigated by grain contact and temperature, 
could contribute to the supposed lower crustal or upper mantle high electrical 
conductivities previously explained ( 28 ) as due to hydrous minerals or partial 
melting. The mafic granulite character, specific gravities, and the T, P and 
Vp estimates for the granulites- correspond to model parameters (25) of the lower 
continental crust, but the extent to which native iron is present requires to 
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be established as does the precise mechanism thpt induced decomposition of 
iron-rich alnandine garnet and il&enite. Decomposition is not due to the 
ldmberlite event because these gas-charged CC^-rich eruptives are moderately 
oxidizing, and because mantle metasomatism is typically potassic (30), whereas 
the native iron-bearing granulite xenoliths are highly reduced and are sodium- 
em iched. 
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Table 1. Magnetic properties and specific gravities of lower crustal granulites. 


Sample 

(a) 

SG (b) 

\ <c> 

(S.I. 

xl0' A ) 

NEM (d) 
(A.nf 1 ) 

J (e) 

(A.nf 1 

xl0“ 2 ) 

<*N 

O' 

M (g) 
(A.nf 1 ) 

1 

3.12 

26 

0.21 

6.5 

3.2 

0.28 

3 

2.75 

7.4 

0.0062 

1.8 

0.34 

0.024 

4 

3.21 

13 

(±3) 

1.6 

3.2 

50 

1.6 

7 

3.02 

8.2 

0.009 

(±0.001) 

2.0 

0.45 

0.029 

8 

3.01 

8.6 

0.19 

2.1 

9.0 

0.21 

9 

3.18 

18 

(±2) 

9.1 

(±0.5) 

4.5 

200 

9.1 

12 

3.06 

15 

0.96 

3.7 

26 

1.0 

13 

2.87 

5.3 

(±1) 

0.073 

(±0.004) 

1.3 

5.6 

0.086 

14 

2.97 

8.4 

0.14 

2.1 

6.7 ■ 

0.16 

17 

3.25 

11 

0.085 

2.7 

3.1 

0.11 

26 

3.08 

5.7 

0.0043 

1.4 

0.31 

0.018 

(a) Series SC-C; 2.54 

x 2.54 cm 

cylindrical 

cores; measured at room 

temperature 


and pressure, (b) SG = apparent specific gravity, (c) X^ = apparent magnetic 
susceptibility per unit volume; error is ±5 x 10”^ except v/here noted, (d) NRM 
= natural remanent magnetization; error is ±2 in L.S.L, except where noted. 

(e) J = induced magnetization, calculated as J = X .H where H . for 
Liberia = 24.9 A.nf 1 . (f) Q = Konigsberger ratio = NEM.J -1 . (g) M a total 
magnetization = MEM + J. 
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Figure Captions 

Fig. 1. Reflected light phot omicrographs of native iron-bearing assemblages 
using a x20 oil-immersion objective: (a) shows garnet decomposition along 

cracks to metallic iron (bright, white spheres) in association with scapolite 
+ tschermakitic amphibolej (b) is an ilmenite grain, set in clinopyroxene, 
that has undergone subsolidus reduction to lamellar ulvospinel (dark gray 
laths) + iron metal (bright, white), with the black laths being plucked areas 
from sample preparation. 

Fig. 2. Relation of natural remanent magnetization (NRK) to apparent magnetic 
susceptibility per unit volume (X y ) for lower crustal granulites (from Table 1). 




